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X-ray (Mo K alpha) studies were made of pure benzene, 
pure cyclohexane, and their mixtures, using a Geiger- 
Mueller counter with balanced filters. The benzene showed 
at 25°C a rather sharp intense peak at the customary 8°40’ 
and a new broader peak of about one-eighth intensity at 
4°20’. The cyclohexane showed a similar pair of peaks at 
8°0' and 4°0’. The mixtures all showed a broad peak of low 
intensity reaching from 4°0’ to 4°20’ and three sharp peaks 
at 8°0’, 8°20’, and 8°40’ superimposed on a broad peak 
which reached from about 8°0’ to about 8°40’. The positions 
of the three sharp peaks were not affected by the pro- 
portions of the two components, and their relative in- 
tensities were independent of amount of stirring. The 


customary shift with temperature was observed on all 
peaks. The results are consistent with a cybotatic type of 
picture. Benzene is pictured as having short-lived highly 
organized groups of -B-B-B-B- floating in a chaotic sea 
of disorganized benzene molecules. The peak at 8°40’ is 
thought to be due to the presence of vibrating groups 
whose extreme positions simulate groups of -B—B-B— 
B-B—B- and —B-B—B-B—B-B—. Cyclohexane would 
contain similar groups of -C-C-C-C-. The mixtures are 
assumed to contain -B-B-B-B-, -C-C-C-C-, and 
-B-C-B-C-. The presence of any considerable amount 
of more complicated groupings is obviously highly 
improbable. 





INTRODUCTION 


HE structure of liquids has been studied 

by many methods, including that of x-ray 
diffraction,'? but their structure is not yet well 
understood. It is the purpose of this paper to 
give the results of an x-ray study of the mutual 
solutions of benzene and cyclohexane. These 
results, differing from those so far reported in 
the literature,*-* lead to a better understanding 


1G. W. Stewart, Chem. Rev. 6, 483 (1929). 

*G. W. Stewart, Rev. Mod. Phys. 2, 116 (1930). For 
reviews of papers and theories see (a) J. T. Randall, The 
Diffraction of X-rays and Electrons by Amorphous Solids, 
Liquids and Gases (John Wiley and Sons, Inc., New York, 
1934); (b) C. Drucker, Physik. Zeits. 29, 273 (1928). 
fi. Krishnamurti, Ind. J. Phys. 3, 331 (1929); 3, 507 

29). 

‘A. W. Meyer, Phys. Rev. 38, 1083 (1931). 

°R. W. G. Wyckoff, Am. J. Sci. 5, 455 (1923). 
wean E. Murray and B. E. Warren, J. Chem. Phys. 7, 141 
IS9), 

’S. Parthasarathy, Phil. Mag. 18, 90-7 (1934). 
*H. K. Ward, J. Chem. Phys. 2, 153 (1934.) 


of the physical state of liquid structure. All 
authors, with the exception of H. K. Ward, 
seem to have obtained only one peak in their 
diffraction curves from liquid solutions. Ward,® 
working with benzene-cyclohexane solutions, 
obtained two peaks corresponding to those given 
by pure benzene and pure cyclohexane. Murray 
and Warren® repeated Ward’s work and were 
unable to obtain two peaks. Since the experi- 
mental results conflicted with each other and 
with the thermodynamic data of Scatchard, 
Wood and Mochel,’ a new x-ray investigation 
was made on the same system with an entirely 
different type of x-ray diffraction apparatus. 
The new data showed not one or two peaks, but 
four. Since the data were definitely reproducible, 
they seemed to warrant further study. 


® G. Scatchard, S. E. Wood and J. M. Mochel, J. Phys. 
Chem. 43, 119-130 (1939). 
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In the present work a Geiger-Mueller counter 
method was used with compensating filters. This 
method appears to be much more sensitive than 
the photographic method and less cumbersome 
than the ionization-chamber method. 


APPARATUS 
A. X-ray radiation 


A (G.E. X-Ray Corp.) molybdenum anode, 
Coolidge type tube was operated at 42 kilovolts 
(r.m.s.) and a tube current of 30 milliamperes. 
The voltage and tube current were maintained 
to within +1.5 percent by use of voltage 
regulators on the transformer primary circuits 
(Raytheon Manufacturing Company). 

To obtain a monochromatic beam of high 
intensity, the method of balanced filters was 
used.!°!! A filter of filter paper saturated with 
Sr(NO3)2 was made such that its x-ray trans- 
mission was equal to a ZrO, filter (Patterson 
Screen Company) for all wave-lengths except the 
Ka doublet of the tube. This was carefully bal- 
anced and checked over the entire x-ray spectrum 
using the (100) plane of NaCl. 


B. Spectrometer and slit system 


A Spencer spectrometer (No. 818), on which 
angles could be read to one minute of arc, was 
mounted in a rigid frame in such a way that the 
movable arm moved in a vertical plane. Soller 
slits" were substituted for the slits and telescope. 
The slits for the collimator and Geiger-Mueller 
counter consisted of 6 parallel slits 20 cm long, 
0.075 cm wide, by 1.0 cm high. The lead foil 
between slits was 0.0075 cm in thickness, the 
spreaders 0.075 cm. 

The geometry of the incident beam was such 
that only three of the six collimating slits could 
be used. The thickness of liquid used in the 
specimen holder was such that all six of the slits 
could be used on the swinging spectrometer arm. 


C. Sample-holder and temperature control 


To eliminate the possibility of any contamina- 
tion the sample-holder was constructed entirely 


10 P, A. Ross, Phys. Rev. 28, 425A (1926). 
11W. C. Pierce, Phys. Rev. 38, 1409 (1931). 
12 W. Soller, Phys. Rev. 24, 158 (1924). 
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Fic. 1. Sample-holder. 


of glass. (See Fig. 1.) The cell windows were of 
Pyrex glass approximately 0.007 cm in thickness. 
The thickness of the liquid in the cell was 1.7 cm. 
This was large enough to insure a high intensity 
of diffracted beam. It was slightly less than the 
optimum thickness given by 


t=1/n, 


where ¢ is the thickness and uy is the linear 
coefficient of absorption. 

The liquid under investigation was brought 
up to temperature, (+0.5°C), in a glass coil 
immersed in a thermostat. The liquid was 
circulated through the coil and the cell, by 
means of a Pyrex paddle, at a linear speed of 
approximately 3 cm per second. The temperature 
difference between the cell and the thermostat 
was therefore negligible. Temperatures between 
5°C and 10°C were obtained by circulating ice 
water through a cooling coil immersed in the 
thermostat. For a temperature of —21°C a dry 
ice-acetone mixture was syphoned through the 
cooling coil, and a eutectic mixture of NaCl and 
water was used in the bath. 


D. Counting circuit 


The Geiger-Mueller counter was argon-oxygen- 
filled. The cathode had a sufficiently large 
opening to allow all x-rays coming from the slits 
to enter the counter. The voltage supply was 
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maintained at 850 volts by a voltage regulator 
of the Street and Johnson type.” 

The amplifier circuit consisted of three stages 
of voltage amplification, using a 57 and finally 
two 56 tubes. The amplified pulse was fed into 
a modified Pickering recording circuit, using a 
single 885 thyratron tube. The final counts were 
recorded by a Cenco counter in the plate circuit 
of the thyratron tube. With this circuit it was 
found necessary to shield the first stage of 
amplification. The time constants for the circuit 
were sufficiently small to allow the 60 cycle of 
the incident beam to be recorded. A small 
loudspeaker was also built in the circuit, powered 
by a 45 tube, so that the quality of sound from 
the speaker gave an indication of the proper 
operation of the G-M counter and of its circuit. 


MATERIALS 


The benzene was Baker’s “C.P. Thiophene 
Free” with the following physical constants: 
B.P. 80°C; M.P.+5.5°C; index of refraction 
1.5009 at 20°C. 

The cyclohexane was obtained from Rohm and 
Haas Co. It was purified by a method somewhat 
similar to that used by other workers.**® The 
material as received was shaken for 12 hr. with 
fuming H2SO, (25 percent SO3), washed with 
water, shaken with alkaline KMnQ, and finally 
shaken with water again. It was then dried 
over metallic sodium and finally distilled at a 
high reflux ratio in a 15-plate column. Further 
purifications were carried out by fractional 
recrystallizations. The physical constants of the 
final material were: B.P. 80.5°C; M.P.+6.4°C; 
and index of refraction 1.4263 at 20°C. 

The solutions, to be examined by diffraction 
methods, were prepared on a weight basis. They 
had weight compositions of (1) 74.5 percent 
benzene-25.5 percent cyclohexane; (2) 51.4 per- 
cent benzene-48.6 percent cyclohexane; (3) 25.4 
percent benzene-74.6 percent cyclohexane. These 
corresponded to mole fraction compositions of 
(1) 0.76 benzene-0.24 cyclohexane; (2) 0.53 
benzene-0.47. cyclohexane; (3) 0.27 benzene- 
0.73 cyclohexane. The exact percentages had no 
special significance, but were those obtained 


‘8 J. C. Street and T. H. Johnson, J. Frank. Inst. 214, 
155 (1932). 
*W.H. Pickering, Rev. Sci. Inst. 9, 180 (1938). 
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after measuring the approximate amounts by 
volume. The index of refraction was measured 
before and after each run and was found to be 
the same within 0.0001. This indicated that the 
composition did not change during the course 
of the experiments. 


EXPERIMENTAL PROCEDURE 


To be certain that the apparatus was in a 
condition to operate uniformly over a _ long 
period, it was in every case allowed to warm up 
for one hour before making any measurements. 
The temperature of the G-M tube was always 
kept within one or two degrees of 25°C since 
these tubes have a marked decrease in efficiency 
with rising temperature. 

The data reported are for the Ka doublet 
lines of molybdenum, 0.710A. This was obtained 
by subtracting the “strontium filter’’ reading 
from the ‘‘zirconium filter’ reading, the differ- 
ence being the Ka contribution to the diffracted 
beam intensity. The background correction due 
to cosmic rays need not be subtracted because 
its average intensity was eliminated by the two 
filter readings. 

A blank run with the empty sample-holder 
gave a very small diffraction scattering with no 
peaks in the region used in this work. The 
counting rates were never higher than 300 per 
min., so that it was possible to detect failure of 
the mechanical counter by an audible method. 
The loudspeaker volume was adjusted so that 
it was not audible to the operator when the 
counter was operating. Any skipping of the 
counter made it possible to hear the loudspeaker 
pulse and thus to detect the failure of the 
counter. 

The ZrO, and Sr(NOs3)2 filters were changed 
after each minute interval, and all readings were 
thrown out for which any counter failure was 
detected. 

Statistics of Geiger-Mueller counters have 
been developed by several authors.'5~'? Hughes’ 
points out that 


(n—x)’=x, 


1% A. E. Ruark and F. E. Brammer, Phys. Rev. 52, 322 
(1937). 

16. Alaoglu and N. M. Smith, Jr., Phys. Rev. 53, 832 
(1938). 

17 A. L. Hughes, Am. Phys. T. 7, 271-292 (1939). 
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Fic. 2. X-ray diffraction curves; 100 percent benzene. 


where is the number of counts in a given time 
interval and x the average number arriving in 
that time interval. The mean deviation is 
therefore equal to the square root of the counts 
recorded. Sufficient counts were taken to insure 
a valid graph over the range studied. In the 
region between 8°00’ and 8°40’ enough counts 
were made so that the differences between 
counts of the peaks and of the adjacent minima 
between the peaks were always greater than 
twice the mean deviation. 

The shapes of the curves and peak positions 
were determined for each liquid sample at 40°C, 
25°C, and in the case of pure benzene, pure 
cyclohexane and the 51.4-percent benzene mix- 
ture, additional curves were taken at 60°C, 10°C 
and as close as possible to the melting points. 
The melting points are +5.5°C, +6.4°C, and 
— 22°C,}8 respectively. 

The shapes of the curves were reproducible 
but the relative intensities varied from week to 
week. This was probably due to long time 
changes in the characteristics of the counter 
circuit. In order to put all curves on a common 
intensity scale, the following method was used: 
(1) Data on all the peaks of all the liquid 


18 Int. Crit. Tab., Vol. IV, p. 133. 
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Fic. 3. X-ray diffraction curves; 100 percent cyclohexane. 


samples studied were taken at each of the 
constant temperatures mentioned above, using 
a separate day for each temperature. (2) Meas- 
urements of the peaks of the 51.4-percent 
benzene mixture (M.P. —22°C) at —21°C, 10°C, 
25°C, 40°C and 60°C, were made on the same 
day. (3) The data from (2) were used to correct 
the measurements of (1) so that the final results 
shown in Figs. 2-5, were what might have been 
expected if all the data could have been obtained 
on a single day. 


EXPERIMENTAL RESULTS 


Diffraction curves for the pure liquids and 
mixtures at various temperatures are shown in 
Figs. 2-5. One unit on the intensity scale 
represents one hundred counts per five minutes. 
26 is the angle of deviation from the incident 
beam made by the diffracted x-rays. 

The range of 26 studied was from 2°00’ to 
11°00’. This is not entirely reproduced in the 
curves shown, because of the size of graphs 
necessary. In all cases a sniall intensity peak 
was observed at a low angle of diffraction. The 
exact position of this peak was carefully checked 
in several cases. In the case of pure benzene 
and cyclohexane it was found, within the 
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Fic. 4. X-ray diffraction curves; 51.4 percent benzene 
solution. 


accuracy of measurement, to be at one-half the 
angle of the high intensity peak (for small angles 
sin @=8@). For the mixtures no sharp peak was 
obtainable but a diffuse peak was found between 
26=40°00 and 20=4°20. The intensities of these 
peaks were approximately one-eighth of the high 
intensity peaks. 

The lowest temperature curves in the cases of 
Figs. 2-4 were taken to determine whether the 
liquid structure might change at a temperature 
just above the melting point. 


INTERPRETATION OF EXPERIMENTAL RESULTS 


Several theories have been proposed to explain 
x-ray diffraction in liquids. The more important 
of these are those due to: Raman and Raman- 
athan,!® Zernike and Prins,2°-2 Debye® and 
Stewart.2, The approach of the first three is 
definitely mathematical in character, while that 
of Stewart is a nonmathematical picture of the 
state of the liquid. Raman and Ramanathan’s 
function was able to predict only one peak in 
—9C, V. Raman and K. R. Ramanathan, Proc. Ind. 
Assn. for Cultivation of Science, p. 127 (1923). 

1927) Zernike and J. A. Prins, Zeits. f. Physik 41, 184 

a A. Prins, Zeits. f. Physik 56, 617 (1929). 


J. A. Prins, Naturwiss. 19, 435 (1931). 
*° P. Debye, Physik. Zeits. 31, 348 (1930). 
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Fic. 5. X-ray diffraction curves; 74.5 percent benzene, 74.6 
percent cyclohexane solutions. 


the diffraction curve and its exact position had 
to be obtained experimentally. Zernike and 
Prins used a distribution function evaluated in 
an empirical way. Debye used the experimental 
data to determine the distribution function. 

These methods were applicable in the case of 
molecules having spherical symmetry, and many 
liquids have been worked out on the basis of the 
theories of Prins and Debye. The application of 
these theories to unsymmetrical molecules such 
as benzene or cyclohexane has not been so 
successful. Several distribution functions have 
been used for benzene.**~*’ All of these functions 
gave fair checks at high diffraction angles, but 
not at angles less than the principal peak. It 
seems that one must resort to the method of 
Stewart to explain the experimental facts in the 
cases of solutions. 

It is well accepted that a liquid approximates 
a crystalline powder in the sense that the 
diffraction peak positions must be close to those 
given by Bragg’s formula 


nd = 2d sin 86. (1) 


*S. Katzoff, J. Chem. Phys. 2, 841 (1934). 

*% L. Bewilogua, Physik. Zeits. 33, 688-692 (1932). 
26 Gustav Thomer, Physik. Zeits. 38, 48-57 (1937). 
27C. M. Sogani, Ind. J. Phys. 1, 357 (1927). 
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TABLE I. Peak positions. 
—-C-C-C- —-B-B-B- —-B-C-B-C- 
COMPOSITION TEMPERATURE 2 d 2 d 2 d 

100% CeHi2 6.7°C 8°03’ 5.064 

4 10°C 8°02’ 5.07 

a = 25°C 8°00’ 5.09 

5 5 40°C 7°58’ 5.11 

38 e 60°C rie 5.14 
100% CeHe S.7°C 8°42’ 4.68A 

" ‘a 10°C 8°41’ 4.69 

3 ni ys gg 8°40’ 4.70 

say 6 40°C 8°38’ 4.72 

= a 60°C 8°35’ 4.74 
51.4% CeHe —21°C 8°04’ 5.05A 8°43’ 4.67 8°24’ 4.85A 

- - 10°C 8°02’ 5.07 8°42’ 4.68 8°22’ 4.87 

= . 25°C 8°00’ 5.09 8°40’ 4.70 8°20’ 4.88 

i 9 40°C 7°58’ 5.11 8°38’ 4.72 8 18’ 4.90 

= 60°C 7°56’ 5.13 8°36’ 4.73 8°15’ 4.93 
74.5% CoHe 25°C 8°00’ 5.09 8°40’ 4.70 8°20’ 4.88 

‘i ws 40°C 7°58’ Sf 8°39’ 4.71 8°19’ 4.89 
74.6% CoeHie Fg 8°00’ 5.09 8°40’ 4.70 8°20’ 4.88 

iz a 40°C 7°59’ 5.10 8°38’ 4.72 8°18’ 4.90 
62.4% CeHe 25°C 8°00’ 5.09 8°40’ 4.70 8°21’ 4.88 
62.8% CeHiz aa 8°00’ 5.09 8°40’ 4.70 8°19’ 4.89 








Applying this formula (1) to the data of Figs. 2 
and 3, we obtain, for the x-ray radiation used 
(\=0.710A), the following values of d at 25°C: 
for 100 percent cyclohexane, 20=8°00' and 
d=5.09+0.01A; for 100 percent benzene, 26 
= 8°40’ and d=4.70+0.01A. Such values of d 
have been interpreted by several authors as the 
effective thickness of the molecules in the 
liquid.”8 

It remains to explain, if possible, the low peaks 
at 26=4°00’—4°20’, and the three peaks shown 
by all the mixtures at 20= 8°00’, 8°20’ and 8°40’. 
Let us consider the solution to be made up 
partly of randomly distributed molecules of 
benzene and cyclohexane, and partly of small 
groups, each existing temporarily, whose mole- 
cules have an almost orderly arrangement. This 
is consistent with Stewart’s explanation’® of 
x-ray diffraction in pure liquids. It is inherent 
in the theory of x-ray diffraction that sharp 
diffraction peaks require the diffracting material 
to have a periodic structure, with enough 
repetition of the fundamental structure to give 
a noticeable contribution to the diffraction 
pattern. For convenience of reference we shall 
use the symbol (—B—B-B-) to represent some 
definite grouping, composed entirely of pure 
benzene, the symbol (-C—C—C-) to represent a 
similar. grouping of pure cyclohexane, and 


28 G. W. Stewart, Phys. Rev. 33, 889 (1929). 


the symbols (-B—C-—B-—C-), (-B—C—-C-—B-C--C-), 
(-B-C—B-C-C-B-), etc., to represent groups 
composed of both benzene and cyclohexane. 
Obviously the groupings (—B-B-B-) and 
(-C-C-C-) show a repetitive structure with the 
single molecule as the fundamental unit. We 
should therefore expect definite diffraction peaks 
corresponding to the grating spaces given by 
adjacent molecules. Such peaks ‘were observed 
experimentally in the present work. Since we 
are dealing with nonrigid liquids, instead of 
solid crystals, we must expect the pure liquid to 
show an equilibrium state between (—B,—B,— 
B,-B,-) and (-B,—B,--B,—B,-) or between 
(-C,-C;—C,,-C,,-) and (-C.—C -C,,—C,,-). This 
should give a diffraction peak corresponding to 
the distance, B,-B,, (or C.-C,,). This spacing is 
obviously not far from twice the mean spacing 
between adjacent molecules. Such peaks have 
been observed experimentally in every case. 
They are not shown in Figs. 2—4 for lack of space. 

The groupings (—B,-C,-B;-C,-) show a 
slightly more complicated periodic structure. 
Such a grouping should show a diffraction peak 
corresponding to the fundamental spacing 
-B,-B,;- or —C,-Cy-, but since the diffracting 
power of B and C are similar, this peak should 
be weak due to interference between B and C. 
Such peaks have been observed in the case of all 
the mixtures studied. The spacing —B,-C.- 
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obviously should be intermediate between the 
spacings —B—B- and —C—C- appearing in the 
(—B-B-B-) and (-C-C-—C-) groupings, and this is 
found to be the actual case. Because of the 
similarity in diffracting power of B and C, a 
definite diffraction peak should be found for all 
mixtures of B and C at a distance corresponding 
to -B,-C.-, (ie., half the spacing B,-B,). 
Experimentally such peaks were found for all 
our mixtures. They are shown in Figs. 4 and 5, 
and in detail in Table I; within experimental 
error, these peaks fall half-way between the 
peaks of (-B—B-B-—) and (—C-C-C-). 

The next most complicated grouping can be 
represented by (—B.-C,;-C,-B,-C-C;-). The 
basic repetitive distance for this grouping is 
BB, This distance would have to be repeated 
several times to give a definite diffraction peak. 
The probability of many successive units of 
(-B-C-C-) occurring in a single group is so 
small that any diffraction peak would necessarily 
be weak. Even if such a definite peak were 
strong enough to be detected, it would lie at too 
small a diffracting angle to be observed on our 
apparatus with Mo Ka rays. A similar argument 
could be used for still more complicated mo- 
lecular groups. 

It might be argued that the shapes of the 
experimental curves are due to addition of two 
complete curves for (-C-C-C-) and (—B-B-B-). 
Figure 6 shows, by dotted lines, the curves for 
(-B-B-B-) and for (—C-C-C-) obtained by 
multiplying the experimental curves (at 25°C) 
of Figs. 2 and 3 by the mole fractions 0.534 and 
0.466, respectively. This composition corresponds 
to the weight fractions 0.514 and 0.486 shown 
in Fig. 4. Figure 6 shows that the sum of the 
two dotted line curves does not even approxi- 
mate the experimental curve. 

Considering the diffraction effects of such 
groups as (—B,-C;-C,-B,-C-C;-) mentioned 
above, we see that scattering over an angular 
range ought, however, to occur due to spacings, 
B-C;, Cr-C,, C,-Br, etc., because of the 
similarity of diffracting power of B and C, and 
the close similarity of the distances, B,—Cy, 
C;--C,, C,-By, etc. The same argument applies 
also to more disordered groups, such as 
(-B-C-—B-C-—C-B-), and to molecules which may 
be in such a chaotic configuration as not to show 
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even a definite molecular orientation. The 
contribution to the diffraction curve due to all 
such groupings should be a diffuse peak, with 
its maximum between 2@=8°00' and 8°40’. 

To summarize the above, we may consider 
that the total x-ray diffraction curve is made up 
of the following:—(1) sharp peaks due to the 
order states (—B-B-B-), (-—C-C-C-), and 
(-B—C-—B-C-) and (2) a diffuse peak, super- 
imposed on the peaks of (1), due to more 
complicated groupings and randomly distributed 
molecules. 

It remains to determine the relative magnitude 
of the order and disorder contributions to the 
diffraction curve, both in the case of the pure 
liquids and in the case of the solutions. This 
may be done as follows:—(1) A smooth curve 
was drawn on the curves of Figs. 2-5, connecting 
the points on both sides of the main peak, 
where the slope seemed to change most rapidly. 
This was entirely arbitrary but is apparently 
justified by the end results. It is similar to the 
method of Stewart in obtaining half-peak 
width.”’ The area under this curve was considered 
to represent the contribution due to all the 
disorder states. (2) The area above the smooth 
curve of (1) was considered to be due to the 
x-ray diffraction of the order states (-B—B-B-), 
(-C-C-C-) for pure benzene and pure cyclo- 
hexane, respectively, and of the order states 
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Fic. 6. Theoretical curve due to addition of curves for 
benzene and cyclohexane. Dotted line curves from pure 
benzene and cyclohexane curves of Figs. 2 and 3. (C+B): 
Addition of dotted line curves. Upper curves: Actual 
experimental curve. 








St 
UIA iaieeestes suid 


ee ene ead 


a= 





448 P. H. BELL AND W. P. DAVEY 


INTENSITY 





2c 


Fic. 7. Theoretical diffraction curves for 51.4 percent 
benzene solution. 


(—B-—B-B-), (-C-C-—C-) and (-B—C-B-C-) of the 
mixtures. If we assume that the (-B—B-B-) and 
(-C-C-C-) groupings in the mixtures are in 
every way the same as in the pure liquids, then 
the amounts of (—B-B-B-) and (—C-C-C-) 
present in a given mixture should be proportional 
to the mole fraction of each in that mixture. 
In other words, the intensity due to (-C-—C-—C-) 
in the pure cyclohexane times the mole fraction 
of C in the solution should equal the intensity 
due to (-—C-—C-C-) in the solution. A similar 
argument should hold for benzene. Using the 
(—B-—B-B-) and (—C-—C-—C-) values of Figs. 2 and 
3 placed on an arbitrary base line, as shown, the 
theoretical curves of Figs. 7 and 8 were obtained. 
The (—B—C—B-C-) peak intensity was made such 
that the addition of all curves would give the 
experimental peaks of Figs. 4 and 5. 

Using the method of (1) and (2) the peak 
intensities were determined and tabulated in 
Table II. The intensities are expressed in terms 
of counts per 5 minutes, as before. Those given 
in parenthesis are the theoretical values from 
Figs. 7 and 8. The others are from Figs. 2-5. 
The half-peak widths are for the order state 
curves for the pure liquids, and over the entire 
three-order peaks for the solutions. 


The complete theoretical curves are shown in 
Figs. 7 and 8. The top dotted line is the addition 
of all the curves contributing to the x-ray 
diffraction. Curves of Figs. 7 and 8 and those of 
Figs. 4 and 5 are in good agreement. The data 
for absolute intensities in the case of the 51.4 
percent benzene mixture are better than for the 
other mixtures since they were very carefully 
checked. The other mixtures were studied more 
to prove whether all the peaks could be detected 
than to ascertain their absolute heights. Phe 
theory did not give very good agreement in the 
case of the —21°C curve of the 51.4 percent 
benzene mixture, since it was necessary to use 
the data for the pure liquids, approximately 15 
degrees higher. 

A few possible objections to this present work 
must be met and, if possible, satisfied. The 
matter of peak resolution, which was first 
considered by Meyer,?® and later by Murray 
and Warren,® seems to be adequately answered, 
if the interpretation given in Figs. 7 and 8 is 
correct. It is true that these separate peaks have 
not been observed by other workers, but most 
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Fic. 8. Theoretical diffraction curves for A, B—74.5 
percent benzene; C, D—74.6 percent cyclohexane solu- 
tions. 


29 A. W. Meyer, Phys. Rev. 38, 1083-1093 (1931). 
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TABLE II. Relative peak intensities. Values in parentheses are from Figs. 7 and 8. All other values are experimental data 
shown in Figs. 2-5. 








COMPOSITION TEMPERATURE 


HALF-PEAK 
WIDTH 


RELATIVE PEAK INTENSITY 





6.7°C 
10°C 
25°C 
40°C 
60°C 
5.7°C 
10°C 
25°C 
40°C 
60°C 
—21°C 
10°C 
25°C 
40°C 
60°C 
25°C 
40°C 
25°C 
40°C 
25°C 
25°C 


100% CcHiz 


74.5% CoH 
74.6% CcHiz 


62.4% CeHe 
62.8% CeHie 


30’ 

35’ 

35’ 

35’ 

25’ 

20’ 

30’ 

30’ 

30’ 

30’ 

(140) 55’ 
(110) 1°05’ 
(90) 1°00’ 

(80) 1°05’ 

(60) 1°05’ 

(100) 55’ 
(90) 1°05’ 

(130) 1°00’ 
(100) 1°00’ 
(130) 55’ 
(110) 1°00’ 


(130) 
(110) 
(85) 
(70) 
(60) 
(110) 
(100) 
(40) 
(35) 
(90) 
(60) 








of them except Ward* used a_ photographic 
method. It seems questionable that a photo- 
graphic method would ever resolve peaks of such 
intensities as those obtained, because of the high 
background density of the film, and the small 
distance on the film between peaks. The experi- 
mental peaks, it is true, as pointed out by 
Murray and Warren, should not be sharp but 
slightly rounded. The data of the present work, 
taken at small angular intervals show a definite 
rounding, although the scale of Figs. 2-5 does 
not bring this out clearly. Also one might ask 
whether the (-B-—B-—B-) and (—C-C-C-) group- 
ings found in the mixtures were due to insuffi- 
cient mixing of the components. Experimentally, 
the x-ray diffraction results were independent of 
the amount of stirring. The effects observed 
could not have been caused by a streaming effect 
in the liquid, for curves run at room temperature, 
with and without stirring, were identical. 

The results reported here may be discussed 
with relation to other data from the literature. 
The data for the unit diffracting distance in a 
given group, (d=n/2 sin 6), read off a graph 
of d versus T, over a temperature range from 
10°C to 60°C (see Table I), show that the 
increase in d is more than 6 percent and less 
than 18 percent, with a probable value of 12 
percent. This is consistent with the results now 


in the literature***°-* on x-ray diffraction of 
liquids. The fact that the I.C.T. tables show, 
for both benzene and cyclohexane, a volume 
expansion of only 6 percent, indicates strongly 
that only a small fraction of the liquid can be 
in an order state at any one time. This conclusion 
is consistent with thermodynamic data on liquids 
in general,**%* and on the benzene-cyclohexane 
system.’ The positions of the (-B—B-B-) and 
(-C—C-C-) peaks at the freezing points of the 
pure liquids strongly indicate that the order 
state is nearly that of the solid, in one dimension 
at least. The two strong lines of the solid benzene 
arising from the (111) and (020) planes give an 
average d of 4.65A,*° the experimental value 
from the liquid peak at 5.7°C gives d=4.68A. 
Two very strong lines for solid cyclohexane gave 
an average d of 5.04A** compared to a value at 
6.7°C of 5.06A for the liquid cyclohexane. 

The positive deviation of the benzene-cyclo- 
hexane system from Raoult’s law® and the 
increase in molal volume on mixing®*? would 


30 E. W. Skinner, Phys. Rev. 36, 1625-30 (1930). 


( 928} S. Ramasubramanyan, Ind. J. Phys. 3, 137-49 
1 , 

# V. I. Vaidynathan, Ind. J. Phys. 5, 501-524 (1930). 

33 J. G. Kirkwood, J. Phys. Chem. 43, 97 (1939). 

%G. Scatchard and W. J. Hamer, J. Am. Chem. Soc. 
57, 1805 (1935). 

3% E. G. Cox, Proc. Roy. Soc. A135, 491 (1932). 

86 K, Lonsdale and H. Smith, Phil. Mag. 28, 614 (1939). 

37 H. Poltz, Zeits. f. physik. Chemie B32, 243-273 (1936). 
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lead to the conclusion that the intermolecular 
forces between benzene and cyclohexane mole- 
cules were less than those between benzene- 
benzene or cyclohexane-cyclohexane molecules. 
Therefore, the probability of the (-B—~C—B-—C-) 
array appearing in the liquid should be less than 
that predicted on the basis of mole fractions. If 
the forces had been all equal, the (-B—C-—B-C-) 
peak intensity should have been between the 
(-B-B-B-) and (—C-C-C-) intensities. In the 
case of the carefully studied 51.4-percent benzene 
mixture, Table II shows that not only is the 
(-B-C-B-C-) peak intensity lower than the 
average intensity of the benzene and cyclohexane 
peaks, but also that it is almost always lower 
than the benzene peak itself. This is consistent 
with the above prediction based on Raoult’s 
law. 


CARL GAMERTSFELDER 


SUMMARY 


(1) A low intensity peak at half the diffraction 
angle of the main peaks has been detected and 
a plausible explanation has been given for its 
existence. 

(2) The main peak from x-ray diffraction 
measurements of benzene-cyclohexane solutions 
has been resolved experimentally into three peaks. 

(3) The observed diffraction peaks in the 
liquid mixtures are explained on an_ order- 
disorder basis, which is consistent with existing 
thermodynamic data. 

(4) Theoretical curves are shown to give 
intensities comparable to measured values, for 
the most probable order states in the solutions. 

(5) The usefulness of Geiger-Mueller counters 
for liquid x-ray diffraction studies has been 
demonstrated. 
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Monochromatic x-ray diffraction patterns of liquid chlorine, tin, indium, zinc, cadmium, 
aluminum, and lithium were obtained, using Mo Ka radiation. The effect of temperature change 
on the patterns was determined for tin and indium. Fourier analyses of these patterns led to the 
determination of atomic distribution curves for each element. Each atom in liquid chlorine has 
one nearest neighbor, which means that chlorine is diatomic in the liquid state. The distribution 
curves for liquid tin, zinc and aluminum agree very well with the distributions in the corre- 
sponding crystals. The distribution curves for indium, cadmium, and lithium do not agree very 
well with distributions in the corresponding crystals. 


INTRODUCTION 


TOMIC distributions in liquids are obtain- 
able from x-ray diffraction patterns by a 
method of analysis based on the theories of 
Zernike and Prins! and of Debye and Menke.’ 
Such determinations have been made for the 
elements mercury,”* gallium,* sodium,**® potas- 


1F. Zernike and J. Prins, Zeits. f. Physik 41, 184 (1927). 
931) Debye and H. Menke, Erg. d. Tech. Rintgenkunde II 
1931). 
3R. N. Boyd and H. R. R. Wakeham, J. Chem. Phys. 
7, 958 (1939). 
4H. Menke, Physik. Zeits. 33, 593 (1932). 
5L. P. Tarasov and B. E. Warren, J. Chem. Phys. 4, 
236 (1936). 
( 938). Trimble and N. S. Gingrich, Phys. Rev. 53, 278 
1 " 


sium,’ phosphorus,® sulfur,? and argon,!®" all in 
the liquid state. In the present work the diffrac- 
tion patterns of liquid chlorine, tin, indium, zinc, 
cadmium, aluminum, and lithium have been ob- 
tained and their atomic distribution curves have 
been determined. The chlorine was liquid under 
its own vapor pressure at room temperature, the 


7C. D. Thomas and N. S. Gingrich, J. Chem. Phys. 6, 
411 (1938). 

®C. D. Thomas and N. S. Gingrich, J. Chem. Phys. 6, 
659 (1938). 

9N.S. Gingrich, J. Chem. Phys. 8, 29 (1940). 

10 A. Eisenstein and N. S. Gingrich, Phys. Rev. 58, 307 
(1940). 

11K. Lark-Horovitz and E. P. Miller, Nature 146, 459 
(1940). 
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Fic. 1. (a) Corrected experimental intensity curve for liquid chlorine at 25°C. (b) Radial 
density distribution of atoms about any one atom in liquid chlorine at 25°C. 


metals were just above their melting points, and 
in the cases of tin and indium, temperatures 
above the melting points were also used. 


EXPERIMENTAL 


Two cylindrical cameras were used in obtaining 
the diffraction patterns, one with a radius of 8.85 
cm and the other with a radius of 9.20 cm. Both 
radii were mechanically measured and they were 
also checked by calculations from the diffraction 
patterns of powdered rocksalt. Mo Ka radiation 
monochromatized by reflection from a rocksalt 
crystal outside the camera was directed through 
a pair of collimating slits before being scattered 
by the sample at the center of the camera. Each 
camera was provided with a sample-holder which 
could be heated while the rest of the camera was 
water-cooled to protect the photographic film 
used to detect the scattered x-rays. The liquid 
chlorine was contained in a thin-walled glass 
capillary, 0.15 cm in diameter, under its own 
vapor pressure of about 7.7 atmospheres at 25°C. 


Tin (at 250°C and 390°C), indium (at 160° and 
390°C), zinc (at 460°C), and cadmium (at 
350°C) were flat samples, about 0.008 cm thick, 
clamped between mica sheets 0.0005 cm thick. 
The aluminum (at about 700°C) in the form of a 
wire 0.17 cm in diameter was heated in a small 
Bunsen flame which kept the aluminum liquid 
inside its very thin but tough oxide coat. The 
lithium sample (at 200°C) was shaped with its 
horizontal cross section as a right isosceles 
triangle, and the x-ray beam was incident perpen- 
dicular to the hypotenuse. The length of path of 
the x-rays in the lithium was 1.8 cm. The lithium 
was contained in a special brass holder a1. 1 was 
prevented from reacting with the brass by means 
of a thin coat of carbon applied with a thin 
shellac binder and then baked. Mica windows 
were used with this holder. 

The intensity of scattered x-rays at various 
angles was determined from the blackening of the 
photographic film by means of a microphotome- 
ter. The intensity in each case was then plotted 
as a function of sin 6/A, where @ is half the angle 
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Fic. 2. (a) Corrected experimental intensity curves for liquid tin at 250°C 
and 390°C. (b) Radial density distributions of atoms about any one atom in 


liquid tin at 250°C and 390°C. 
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Fic. 3. (a) Corrected experimental intensity curves for liquid indium at 
160°C and 390°C. (b) Radial density distributions of atoms about any one 
atom in liquid indium at 160°C and 390°C. 


of scattering and \ is 0.710A. The experimental 
intensity curves have been corrected and are 
given as the (a) parts of Figs. 1 to 7. The nature 
of the corrections is given in the next section. 
Table I gives the positions, in sin @/X units, of 
the intensity maxima of the various patterns. 


Prins,!* and Danilov and Radchenko" have ob- 
tained the first maximum in the intensity pattern 
of tin and their results agree very well with the 


2 J. Prins, Trans. Faraday Soc. 33, 110 (1937). | 
13 J. V. Danilov and V. : Radchenko, Physik. Zeits. der 
Sowjetunion 12, 313 (1937). 
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Fic. 4. (a) Corrected experimental intensity curve for liquid zinc at 460°C. 
(b) Radial density distribution of atoms about any one atom in liquid zinc at 


460°C. 
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Fic. 5. (a) Corrected experimental intensity curve for liquid cadmium at 
350°C. (b) Radial density distribution of atoms about any atom in liquid 


cadmium at 350°C. 


ones given here, while Sauerwalde and Teske™ 
report a value which is much larger. Randall!® 


* F. Sauerwalde and W. Teske, Zeits. f. anorg. Chemie 
210, 247 (1933). 

* F. R. Randall, The Diffraction of X-Rays and Electrons 
by Amorphous Solids, Liquids and Gases (John Wiley and 
Sons, Inc., New York, 1934), p. 134. 


reports that the diffraction maximum of liquid 
aluminum corresponds to that of the (111) planes 
of the solid, but no exact data are given. The 
spacing of the (111) planes in aluminum is 
2.33A and this corresponds to a diffraction 
maximum, using Bragg’s law, at sin 6/A=0.214 
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Fic. 6. (a) Corrected experimental intensity curve for liquid aluminum at 700°C. (b) Ra- 
dial density distribution of atoms about any one atom in liquid aluminum at 700°C. 
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Fic. 7. (a) Corrected experimental intensity curve for liquid lithium at 200°C. 
(b) “ioe density distribution of atoms about any one atom in liquid lithium at 
200°C. 
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which is in good agreement with the main peak 
reported here. 

In the course of this investigation it was 
necessary to determine the densities of liquid 
indium and liquid lithium. The density of indium 
at various temperatures was determined by 
filling a bulb on the end of a capillary tube with 
indium. The volume of the indium was measured 
by calibrating the bulb and tube with mercury 
the density of which is very well known. The 
density of liquid indium at its melting point of 
155°C is 7.03+0.01 g/cm* and the coefficient of 
volume expansion is 0.00013 per degree C. The 
density of liquid lithium was found by means of a 
pyknometer filled with oil in which the lithium 
was placed. The density of liquid lithium at 
195°C is 0.46+0.01 g/cm’. 


THE FourRIER METHOD 


The theory, upon which the analyses of these 
diffraction patterns are based, is given elsewhere,'® 
and the equation used here is 


2R 
4 R?p(R) = 42 R2p9+— f Si(S) sin RSdS, (1) 
Tv 0 


where R is the distance in angstroms from any 
arbitrary atom taken as the origin. p(R) is the 
actual density in atoms per cubic angstrom, po is 
average density in atoms per cubic angstrom, and 
S=4n sin 6/X. The function i(S) is defined as 
(1/N— F*)/F? where I/N is the intensity of the 
coherent scattered radiation per atom from a 
sample having no absorption and with correc- 
tions for polarization of the scattered x-rays, and 
F is the structure factor of the atom. The values 
of 4rR’p(R)dR give, then, the number of atoms 


TABLE I. 








FIFTH 
PEAK 


THIRD FourRTH 


PEAK 


SECOND 
PEAK 


0.285 y } 0.8 
335 ‘ 
330 
335 
335 
.228 
360 
.210 
.340 


FIRST 
PEAK 


0.122 
174 
175 
179 
177 





BE 

Sn 250° 
Sn 390° 
In 160° 
In 390° 
Zn 460° 115 
Cd 350° .203 
Al 700° 10 

Li 200° .195 


PLT III 








© B. E. Warren, J. App. Phys. 8, 645 (1937). 


TABLE II. 








sin6/d F==fn 
0.0 3.00 
2.20 
Be 
1.47 
1.24 

1.007 

.803 

.634 

.510 











in a spherical shell of thickness dR and radius R 
about afy one atom of the liquid. The experi- 
mentally determined diffraction patterns were 
corrected for polarization’’ and absorption. Ab- 
sorption corrections for samples with flat sides 
were easily calculated and corrections for cylin- 
drical samples are given by Blake.'* Correction 
was also made for the incoherent radiation by 
means of the F and =f,” tables given by Compton 
and Allison’® or interpolations from their values, 
except in the case of lithium whose values were 
not given accurately enough. The F values for 
lithium were calculated from the electron distri- 
bution as given by Fock and Petrashen.®° Table 
II gives the values of the F= >f, and =f,” terms as 
calculated for various sin 6/) values. 

The ratio between the coherent and incoherent 
radiation was determined at a large angle of 
scattering where no effects due to interferences 
appear and then the complete incoherent curve 
was drawn in and subtracted from the intensity 
curve. The F? curve was then drawn in and the 
4(S) values calculated. The values of the integral 
in Eq. (1) were determined by graphical inte- 
gration and the values of 47R’p(R) thus deter- 
mined were plotted against R. 

The results for the various elements are shown 
in the (b) parts of Figs. 1 to 7. The positions and 
numbers of the atoms in the crystal lattices are 
given by the straight lines located at various 
abscissae of the graphs, the height of each line 
representing the number of atoms at that par- 
ticular distance. The dotted lines with the tin 
curve give the atoms for gray tin, the solid 


17S. Katzoff, J. Chem, Phys. 2, 841 (1934). 

18 F.C. Blake, Rev. Mod. Phys. 5, 169 (1933). 

19 A. H. Compton and S. K. Allison, X-Rays in Theory 
and Experiment (D. Van Nostrand, 1935), p. 139 and 781. 

20 V. Fock and M. J. Petrashen, Physik. Zeits. der Sowjet- 
union 8, 547 (1936). 
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TABLE III. 
LIQUID SOLID DENSITY SOLID 

ELEMENT Dist. No. Dist No. Atoms/A?3 Liquip (20°C) 

ci 25° 2.01 0.97 — — 0.0238 1.40 g/cm* — 

Sn 250° 3.38 10.0 3.02 4) 
3.16 2 } white tin 0.0356 6.97 7.31 g/cm? 
3.75 4 
iso Sent 

Sn 390° 3.36 8.9 — -- 0.0350 6.86 — 

In 160° 3.30 8.5 3.24 4 0.0371 7.03 7.28 
3.36 8 

In 390° 3.36 8.4 — — 0.0361 6.84 — 

Zn 460° 2.94 10.8 2.66 6 0.0604 6.51 7.14 
2.94 6 

Cd 350° 3.06 8.3 2.98 6 0.0430 7.98 8.65 
3.34 6 

Al 700° 2.96 10.6 2.86 12 0.0530 2.36 2.70 

Li 200° 3.24 9.8 3.03 8 0.0402 0.46 0.534 
3.50 6 








lines for white tin. Table III gives the positions 
of the first maxima in the distribution curves for 
the liquids and the number of atoms at these 
distances as determined by drawing in sym- 
metrical peaks and measuring the areas under 
them. In addition some data are given on the 
crystalline forms and the densities of the liquids 
and solids. 


DISCUSSION OF RESULTS 


The diffraction pattern of liquid chlorine has 
two distinct maxima and three more fluctuations 
about the F? curve. The distribution curve shows 
definitely that the chlorine is diatomic in the 
liquid state and the interatomic distance of 
2.01A agrees very well with the value of 2.009A 
given by Pauling and Brockway”!.as determined 
from electron diffraction patterns of the gas. The 
rest of the distribution curve is not easy to 
interpret quantitatively, though it seems that on 
the average, one atom of one of the neighboring 
molecules is approximately 3A from the reference 
atom while the rest of the neighboring molecules 
are located at a distance of about 4A. The 
number of neighboring molecules at these dis- 
tances is difficult to determine. The values of the 
distribution function at small values of R are 
very sensitive to the shape of the diffraction 
pattern at large angles of scattering, which is the 
most difficult to measure accurately. This effect 
was noticed in the case of chlorine where an 


L. Pauling and L. O. Brockway, J. Chem. Phys. 2, 
867 (1934). 


incomplete diffraction pattern, without the last 
two maxima, was analyzed. The results in this 
case gave a distribution curve very much like the 
one reported here except in the vicinity of the 
peak at 2.01A. The area under the peak was the 
same but it was broader and moved in to about 
1.8A. In the case of argon® the introduction of 
a fourth peak in the intensity curve did not cause 
any appreciable effect in the first maximum of 
the distribution curve. The effect in this case is 
not the same however, because the first maximum 
for argon is at a larger value (about 3.9A) of R 
than in the case of chlorine. 

Both diffraction patterns of tin have three 
maxima. In the 390°C case the maxima are not so 
high but are broader than in the 250°C case. The 
distribution curve for liquid tin shows that there 
are 10.0 nearest neighbors at 3.38A, which is in 
excellent agreement with crystalline white tin 
where there are 10 nearest neighbors at an aver- 
age distance of 3.34A. There is no such correla- 
tion with the gray tin structure. The distribution 
curve at the high temperature does not have as 
definite maxima as at the lower temperature 
because of the increased heat motion of the 
atoms. The distance to the nearest neighbors is 
not changed but the number is decreased. The 
atoms which leave the first concentration join 
with the secondary concentration which moves in 
to a shorter distance. 

The diffraction patterns of liquid indium at 
160°C and 390°C are very much alike except for 


2N. S. Gingrich, Phys. Rev. 59, 290 (1941). 
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the main peak which is narrower and higher for 
the 160°C case. The third maximum is also much 
stronger for the lower temperature. Liquid indium 
at 160°C has, on the average, 8.5 nearest neigh- 
bors at 3.30A, which is quite different from the 4 
atoms at 3.24A and 8 atoms at 3.36A or 12 atoms 
at an average distance of 3.32A in the crystal. 
The distribution curve does not change much at 
the higher temperature except for the appearance 
of a slight secondary maximum which did not 
show at the low temperature, presumably be- 
cause it was too close to the first concentration. 

Zinc and cadmium have very similar hexagonal 
crystal structures, but the diffraction patterns of 
the liquids are not very much alike. The diffrac- 
tion pattern of zinc is very unusual in that a 
small maximum occurs at a value of sin 6/A which 
is much smaller than that for the main peak. The 
pattern for cadmium is very much like that of 
indium except that the first maximum is at a 
larger value of sin @/A. Liquid zinc has 10.8 
nearest neighbors at 2.94A as compared with 6 
atoms’ at 2.66A and 6 atoms at 2.94A or 12 
atoms at an average of 2.80A in the crystal. 
Cadmium has only 8.3 nearest neighbors at 
3.06A as compared with 6 atoms at 2.98A and 6 
atoms at 3.34A or 12 atoms at an average of 
3.16A in the crystal. 

The diffraction pattern of aluminum is like 
that of zinc in that it has a small maximum 
inside the main peak, so that there are four 
maxima altogether. The distribution curve for 
aluminum shows that there are 10.6 atoms at 
2.96A while in the crystal there are 12 atoms 
at 2.86A. 

The diffraction pattern of liquid lithium is not 
known as accurately as for the other elements 
because the sample used was not small compared 
to the radius of the camera used to obtain it. 
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Another source of error might be in the very large 
incoherent correction for which the recoil factor!® 
could not be neglected. The analysis shows that 
liquid lithium has 9.8 nearest neighbors at 3.24A 
as compared with 8 atoms at 3.03A and 6 atoms 
at 3.50A in the crystal. 

Prins* has attempted to express the atomic 
distribution in a liquid by making a comparison 
with the corresponding crystal, where the atomic 
distribution in the liquid is produced by blurring 
the crystalline distribution. Prins has calculated 
diffraction patterns for atomic distributions ob- 
tained in this way which give good agreement 
with experimentally observed patterns. The 
agreement between the liquid distributions re- 
ported here and the corresponding crystal struc- 
tures is good in the cases of tin, zinc, and 
aluminum. In the cases of indium, cadmium, and 
lithium the agreement is not at all good. Coulson 
and Rushbrooke™ have suggested that there are 
likely to be greater structural changes on melting 
in cases where the volume change is large. How- 
ever, zinc and aluminum both have larger volume 
changes than cadmium, and it is cadmium which 
has the largest structural change. Lithium also 
has a very large volume change on melting but 
the number of nearest neighbors is increased 
rather than decreased. 

The author wishes to thank Professor N. S. 
Gingrich, under whose guidance this work was 
done, for his very helpful criticism and advice 
and for making available certain equipment 
obtained through grants in aid of research from 
the American Academy of Arts and Science and 
from the American Association for the Advance- 
ment of Science. 


23 J. Prins, Physica [3], 3, 147 (1936). 
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A high dispersion study of the OH bands of hydrogen peroxide in the photographic infra-red is 
described. The rotational structure of the two hybrid bands at 10,283.68 cm and 10,291.08 
cm~! is analyzed. The small moment of inertia of hydrogen peroxide is found to be 2.786 Xx 10~*° 
g cm? while the harmonic mean of the two larger moments of inertia is found to be approximately 
33.9 X 10-*° g cm?. The planar cis-configuration is ruled out spectroscopically and it is concluded 
that hydrogen peroxide exists in a non-planar form. It is tentatively decided that the vibrational 
levels of hydrogen peroxide are doubled and that the double minimum potential associated with 
the internal rotation of the OH groups has a high and a low barrier. 





INTRODUCTION 


HE position of the OH groups with respect 
to their relative rotation about the O—O 
bond in hydrogen peroxide is of some interest. 
The OH groups can be imagined as preferring to 
assume the planar cis- or trans-configurations or 
some intermediate non-planar configuration. In 
the third case the molecule will have two 
equivalent enantiomorphic forms and the hinder- 
ing potential for internal rotation will have two 
minima of equal depth and, in general, two bar- 
riers of unequal height. The OH groups undergo 
free rotation, hindered rotation or torsional 
oscillation depending on the heights of the 
barriers of the hindering potential. 

Penney and Sutherland! have proposed a non- 
planar structure for hydrogen peroxide where the 
azimuthal angle? is 100°. Both barriers are given 
as being fairly high. These deductions are based 
on theoretical considerations and are shown to 
give a calculated value of the electric moment 
which agrees well with that obtained from 
measurements of the dielectric constants of dilute 
solutions of hydrogen peroxide in dioxane.® 

Evidence of hindered rotation or torsional 
oscillation should appear in the rotational fine 
structure of the infra-red bands. In the hope of 
obtaining such evidence we have recently studied 


* Arthur A. Noyes Fellow in Chemistry. 
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1W. G. Penney and G. B. B. M. Sutherland, J. Chem. 
Phys. 2, 492 (1934). 

? The angle between the OH groups as seen projected on 
a plane normal to the O—O bond. 
(1932) P. Linton and O. Maass, Can. J. Research 7, 81 


the absorption spectrum of hydrogen peroxide 
vapor in the photographic infra-red. The result 
has been the discovery of an unusual pair of 
bands at 9720 which give some information 
regarding the structure and torsional motion of 
hydrogen peroxide. 


EXPERIMENTAL 


The hydrogen peroxide used in these experi- 
ments had a strength which ranged from 95 to 98 
percent. It was prepared from the commercial 30 
percent solution by distillation under reduced 
pressure. The absorption cell was a six-meter 
Pyrex tube (22 mm diam.) which was heated to 
100°C and through which hydrogen peroxide was 
continuously distilled at a pressure of 90 mm. 
Hydrogen peroxide which had been carefully 
purified was found to decompose very little on 
boiling. It was important to keep the all-Pyrex 
system very clean. Contact of the vapor with 
dirt or even a rough spot on the glass was found 
to cause violent decomposition. 

When a tungsten lamp was used as a light 
source, the exposure time required for the high 
dispersion spectrogram was found to be exces- 
sive, consequently a carbon arc was employed 
as a light source. It was found necessary to run 
the carbon arc in an atmosphere of carbon 
dioxide to eliminate emission lines of the red CN 
bands. This light source reduced the required 
time of exposure by a factor of four; thus the 
spectroscopic plate was amply exposed in eight 
hours. During this time approximately two liters 
of the concentrated hydrogen peroxide were dis- 
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Fic. 1. Microphotometer trace of hydrogen peroxide bands at \9720. 


tilled over. The decomposition was only about 
2 percent. 

Exploration spectrograms were taken on 
Eastman 1Z plates with a glass Littrow spec- 
trograph which covered the region between 7500 
and 12,000A. Only one band at 9720A was found 
and this was then photographed in the first order 
of a twenty-one-foot concave grating spectro- 
graph using Eastman 144Q plates. Measurements 
of the high dispersion plates were made with 
reference to third-order iron arc lines. 


Rotational structure of band 


An examination of the high dispersion spec- 
trogram shows that a hybrid band is present 
which has the Q branches characteristic of both 
the parallel and the perpendicular bands of a 
symmetrical rotator molecule. The appearance 
of this band is analogous to that of the 10,477A 
band of hydrazoic acid.‘ In- the case of hydrogen 
peroxide, as in the case of hydrazoic acid, there 
is a strong dependence of the small moment of 
inertia upon the vibrational state, thus the 
parallel sub-bands are separated and an unusu- 
ally marked convergence is present in the per- 
pendicular component. A second hybrid band 
occurring at a slightly lower frequency also 
appears. This band is similar to the first but is a 
little less intense. The °Qx branches® of the 


*E. H. Eyster, J. Chem. Phys. 8, 135 (1940). 

5In this nomenclature the large letter describes the 
behavior of J in transition, the left-hand superscript gives 
the behavior of K, and the right-hand subscript gives the 
value of K in the lower state. 


parallel component and the ?Qx and *¥Qx 
branches of the perpendicular component of the 
two bands have been identified. The branches of 
the lower frequency band which occur at a dis- 
placement of 7.40 cm have been given primed 
symbols. 

A careful scrutiny reveals that at least a part 
of the band system appears to have a quadruplet 
character rather than a doublet character. Thus 
in the instances where K”’ =3, 4, 5, 6a satellite of 
PQx’ occurs at a shift of about minus 3.8 cm~! 
and also a satellite of ?Qx occurs at a shift of 
about plus 3.8 cm~. These attendant lines are 
weak and somewhat difficult to measure; never- 


TABLE 1. Frequency and identification of lines in the absorp- 
tion spectrum of hydrogen peroxide at 9720. 








FREQUENCY 
cm™! 


FREQUENCY 


IDENTIFICATION cm"! IDENTIFICATION 


PO, 





10,274.76 
276.95 
278.55 
281.39 
282.12 
286.06 
288.82 
290.40 
339.21 
346.51 
351.86 
359.56 
364.44 
371.71 
382.27 
391.82 


10,151.49 
164.69 
168.21 
175.94 
179.65 
188.24 
191.94 
199.25 
203.23 
210.28 
213.96 
221.49 
225.51 
231.10 
235.06 
242.41 
263.44 
269.37 
270.72 
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A high dispersion study of the OH bands of hydrogen peroxide in the photographic infra-red is 
described. The rotational structure of the two hybrid bands at 10,283.68 cm and 10,291.08 
cm~! is analyzed. The small moment of inertia of hydrogen peroxide is found to be 2.785 x 10° 
g cm? while the harmonic mean of the two larger moments of inertia is found to be approximately 
33.9 10-*° g cm?. The planar cis-configuration is ruled out spectroscopically and it is concluded 
that hydrogen peroxide exists in a non-planar form. It is tentatively decided that the vibrational 
levels of hydrogen peroxide are doubled and that the double minimum potential associated with 
the internal rotation of the OH groups has a high and a low barrier. 





INTRODUCTION 


HE position of the OH groups with respect 
to their relative rotation about the O—O 
bond in hydrogen peroxide is of some interest. 
The OH groups can be imagined as preferring to 
assume the planar cis- or trans-configurations or 
some intermediate non-planar configuration. In 
the third case the molecule will have two 
equivalent enantiomorphic forms and the hinder- 
ing potential for internal rotation will have two 
minima of equal depth and, in general, two bar- 
riers of unequal height. The OH groups undergo 
free rotation, hindered rotation or torsional 
oscillation depending on the heights of the 
barriers of the hindering potential. 

Penney and Sutherland! have proposed a non- 
planar structure for hydrogen peroxide where the 
azimuthal angle? is 100°. Both barriers are given 
as being fairly high. These deductions are based 
on theoretical considerations and are shown to 
give a calculated value of the electric moment 
which agrees well with that obtained from 
measurements of the dielectric constants of dilute 
solutions of hydrogen peroxide in dioxane.® 

Evidence of hindered rotation or torsional 
oscillation should appear in the rotational fine 
structure of the infra-red bands. In the hope of 
obtaining such evidence we have recently studied 
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the absorption spectrum of hydrogen peroxide 
vapor in the photographic infra-red. The result 
has been the discovery of an unusual pair of 
bands at 9720 which give some information 
regarding the structure and torsional motion of 
hydrogen peroxide. 


EXPERIMENTAL 


The hydrogen peroxide used in these experi- 
ments had a strength which ranged from 95 to 98 
percent. It was prepared from the commercial 30 
percent solution by distillation under reduced 
pressure. The absorption cell was a six-meter 
Pyrex tube (22 mm diam.) which was heated to 
100°C and through which hydrogen peroxide was 
continuously distilled at a pressure of 90 mm. 
Hydrogen peroxide which had been carefully 
purified was found to decompose very little on 
boiling. It was important to keep the all-Pyrex 
system very clean. Contact of the vapor with 
dirt or even a rough spot on the glass was found 
to cause violent decomposition. 

When a tungsten lamp was used as a light 
source, the exposure time required for the high 
dispersion spectrogram was found to be exces- 
sive, consequently a carbon arc was employed 
as a light source. It was found necessary to run 
the carbon arc in an atmosphere of carbon 
dioxide to eliminate emission lines of the red CN 
bands. This light source reduced the required 
time of exposure by a factor of four; thus the 
spectroscopic plate was amply exposed in eight 
hours. During this time approximately two liters 
of the concentrated hydrogen peroxide were dis- 
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tilled over. The decomposition was only about 
2 percent. 

Exploration spectrograms were taken on 
Eastman 1Z plates with a glass Littrow spec- 
trograph which covered the region between 7500 
and 12,000A. Only one band at 9720A was found 
and this was then photographed in the first order 


of a twenty-one-foot concave grating spectro- 
graph using Eastman 1440 plates. Measurements 
of the high dispersion plates were made with 
reference to third-order iron arc lines. 


Rotational structure of band 

An examination of the high dispersion spec- 
trogram shows that a hybrid band is present 
which has the Q branches characteristic of both 
the parallel and the perpendicular bands of a 
symmetrical rotator molecule. The appearance 
of this band is analogous to that of the 10,477A 
band of hydrazoic acid.‘ In-the case of hydrogen 
peroxide, as in the case of hydrazoic acid, there 
is a strong dependence of the small moment of 
inertia upon the vibrational state, thus the 
parallel sub-bands are separated and an unusu- 
ally marked convergence is present in the per- 
pendicular component. A second hybrid band 
occurring at a slightly lower frequency also 
appears. This band is similar to the first but is a 
little less intense. The °Qx branches® of the 


*E. H. Eyster, J. Chem. Phys. 8, 135 (1940). 

5In this nomenclature the large letter describes the 
behavior of J in transition, the left-hand superscript gives 
the behavior of K, and the right-hand subscript gives the 
value of K in the lower state. 


parallel component and the ?Qx and *¥Qx 
branches of the perpendicular component of the 
two bands have been identified. The branches of 
the lower frequency band which occur at a dis- 
placement of 7.40 cm have been given primed 
symbols. 

A careful scrutiny reveals that at least a part 
of the band system appears to have a quadruplet 
character rather than a doublet character. Thus 
in the instances where K”’ = 3, 4, 5, 6a satellite of 
POs’ occurs at a shift of about minus 3.8 cm~! 
and also a satellite of ?Qx occurs at a shift of 
about plus 3.8 cm~'. These attendant lines are 
weak and somewhat difficult to measure; never- 


TABLE 1. Frequency and identification of lines in the absorp- 
tion spectrum of hydrogen peroxide at 9720. 








FREQUENCY 
cm™! 


FREQUENCY 
cm"! 


10,151.49 
164.69 
168.21 
175.94 
179.65 
188.24 
191.94 
199.25 
203.23 
210.28 
213.96 
221.49 
225.51 
231.10 
235.06 
242.41 
263.44 
269.37 
270.72 
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10,274.76 
276.95 
278.55 
281.39 
282.12 
286.06 
288.82 
290.40 
339.21 
346.51 
351.86 
359.56 
364.44 
371.71 
382.27 
391.82 
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theless, they do seem to be real. The satellites 
have been named ax and bx, respectively. Cor- 
responding branches are not seen on the high 
frequency side of the band centers. 

Figure 1 gives a reproduction of the micro- 
photometer trace of the grating spectrogram 
while Table I lists the frequencies*® and identi- 
fications of all lines which are assigned as Q 
branches and also of the satellites mentioned 
above. 

It is to be noted that the ?Qe, ’Qi, “Qo and ”Q, 
branches are not identified. There is no clear 
evidence of these branches although they are the 
strongest branches in the perpendicular band of 
the symmetrical rotator. This is probably due to 
the fact that the hydrogen peroxide molecule is 
a slightly asymmetric top and consequently the 
Q branches involving the levels where K=0 or 
1 are considerably broadened. The °Q; branch 
involving transitions between levels where K =1, 
however, is observed. One, indeed, might expect 
the perpendicular band Q branches to be smeared 
out before the corresponding parallel band Q 
branches because in the limiting case of the 
symmetric rotator’ we know that the transitions 
involving the lowest J levels contribute most of 
the intensity to the latter branches while the 
higher levels contribute most of the intensity to 
the perpendicular band Q branches. When 
transitions between slightly perturbed levels of 
low J value are the principal contributors to the 
intensity of a Q branch then the branch is sharp 
and well defined. On the other hand when transi- 
tions between strongly perturbed levels of higher 
J value substantially contribute to a Q branch 
then the branch becomes greatly broadened and 
loses its identity. 

The observed spectrum is characterized by 
multifold overlapping of sub-bands so that the 
fine structure is considerably confused and it is 
impossible to analyze the J structure of the 
sub-bands. Regular fine spacing is observed only 
on the high frequency side of the band centers. 
The fine spacing shows convergence towards high 
frequencies, and it falls off from 1.65 cm™ to 
1.59 cm. 

6 The frequencies are based on two independent sets of 
measurements on each of two high dispersion spectrograms. 

7 See the expressions for the intensities of the branches 


of the symmetric rotator bands given by D. M. Dennison, 
Rev. Mod. Phys. 3, 313 (1931). 





Analysis of band structure 


The rotational energy of a symmetric rotator 
may be written: 


F(J, K,v) =D" Xi (IVJ+ IJK) (1) 


The constant term Xo” is omitted in the sum- 
mation, and v represents the set of quantum 
numbers required to specify the vibrational 
state. The first coefficients are related to the 
effective moments of inertia as follows: 


Xo" =A,—3(B.+C,), X10" =3(B.+C,), (2) 
A,=h/8r'cl,’, B,=h/8n'cIp’, 
C,.=h/8rcIc’, y - >B,=C,. 


Hydrogen peroxide seems to be a slightly 
asymmetric rotator with B, almost equal to C, 
so that Eq. (1) very closely represents the 
characteristic energy values of the rotating 
molecule. 

The constants Xo:° and Xo2° are evaluated by 
employing the following relations: 


AF’ (K) 
———- = Xo1°+2X 029M +--+; M=(K+3)?+3 
2K+1 
and 
A2F’’(K) 
— ee ++; M=K?*+1, 


where the term differences are given by: 
AF" (K) =°Qx—?Qx41="Qx— °Qx41 
A2F’'(K) =*®Qx-1—?Qx+1. 
The frequencies of the assigned lines were 
found to be compatible with the above relations 


within experimental error. The best values for 
the constants of the equations were found to be: 


Xo" a 9.23, cm—!, X 02° a << 0.0008; cm~!, 


Similarly the rotational constants for the 
excited state, which we will denote by v=3, were 
found to be: 


Xo1°=8.663cm—!, Xo22= —0.00085 cm. 


An identical treatment of the frequencies of 
the Q branches which are denoted by primed 
symbols gave rotational constants which, within 
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experimental error, were exactly the same as the 
above. The two hybrid bands thus seem to have 
precisely the same rotational structure. 

Since it was impossible to analyze the J 
structure of the bands we were obliged to obtain 
an estimate of the rotational constant X10° 
(which must be very nearly equal to X19*) from 
the regular fine spacing which is clearly resolved 
only on the high frequency side of the band 
centers. We took 1.65 cm~ as the best average 
value of the fine spacing. This is related to the 
rotational constant as follows: 


2X 10° = 1.65 cm 
thus 
X 10° = 0.825 cm—!, 


We believe that the observed fine spacing 
arises from R branches of the parallel sub-bands. 
This being the case the convergence is such that 
the above estimated value of X19° may be a 
little low. 

The center of the band denoted by unprimed 
Q branches is given by: 


G(3) — G(0) - °0x+K?(Xo1°—X 01°) 
+K4(X 02° — X 02°) + si 


The best value for this was found to be 10,291.08 
cm~, Similarly the center of the band denoted 
by primed Q branches was found to be 10,283.68 
cm7}, 

The rotational constants X o;°=9.23; cm and 
X 19° ~ 0.82; cm—! correspond, respectively, to the 
values of I4°=2.78,X10-*° g cm? for the small 
moment of inertia and 2J,°I¢°/(Ip°+Ic°) 
~=33.9X10-° g cm? for the harmonic mean of 
the two larger moments of inertia of hydrogen 
peroxide. 


Structure of hydrogen peroxide 


The harmonic mean of the two larger moments 
of inertia is principally dependent on the length 
of the O—O bond in hydrogen peroxide. Calcu- 
lation shuws that the value given above cor- 
responds to an O—O distance of about 1.48A. 

The fact that the observed bands are of the 
hybrid type shows that the OH groups of 
hydrogen peroxide are not in the planar cis-con- 
figuration. If these groups were so situated, the 
molecule would have a point group symmetry 
Cx, and one would only observe pure perpen- 
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dicular or parallel type bands. Incidentally, we 
also know that the OH stretching vibration 
involved in the observed bands is antisym- 
metric. A symmetric stretching frequency would 
give a pure perpendicular band. 

As shown by Linton and Maass the hydrogen 
peroxide molecule has a fairly large electric 
dipole moment. This rules out the trans-con- 
figuration and, thus, since the cis-form is ruled 
out spectroscopically, the hydrogen peroxide 
molecule must have a non-planar structure. 

We are unable to give an estimate of the 
azimuthal angle? of the OH groups from our 
spectroscopic data although it appears that this 
angle cannot be between 85° and 95°. We con- 
clude this because the absence of the Qs, ?Q,, 
RQ, and *#Q, branches indicates the molecule is 
actually an asymmetric top while a calculation of 
the moments of inertia shows that a hydrogen 
peroxide molecule with an azimuthal angle of 
around 90° is not significantly asymmetric. 

It is interesting to note that Lu, Hughes and 
Giguére® in their recent x-ray studies of the 
“hyperol’’ crystal in these laboratories have 
found the hydrogen peroxide molecule to have 
an azimuthal angle of about 106°. 


Discussion of the band structure 


The bands which have been described and 
analyzed correspond closely to the type of band 
which is characteristic of a spindle-shaped sym- 
metrical top molecule that is fairly rigid. The 
complete rotation of the OH groups about the 
O-—O bond in hydrogen peroxide must be 
strongly hindered. A very different band struc- 
ture would appear if there were a free or a weakly 
hindered internal rotation. In this case the band 
would have a more extended structure, a larger 
Q spacing and, in the event of a weakly hindered 
rotation, an irregular Q spacing. 

The occurrence of two hybrid bands separated 
by 7.40 cm~! must be explained. A visual exam- 
ination of the spectrogram and the micropho- 
tometer trace shows that the low frequency band 
must be more than one-half as intense as the 
other band. The ?Q, and ”Q,’ branches which 


8 Chia-Si Lu, E. W. Hughes, and P. A. Giguére, “Crystal 
structure of the urea-hydrogen peroxide addition com- 
pound, CO(NH:2)-H:202,” to be published in J. Am. Chem. 
Soc. 













































462 7s 2 
occur where the background absorption is weak 
seem to have just about the same intensity. 
Assuming relative intensities of the Q branches 
are given by the matrix elements of the electric 
moment for the symmetric rotator given by 
Dennison,’ we calculate that the perpendicular 
component of the low frequency band is about 
two-thirds as intense as the same component of 
the other band. 

Three possible explanations for the hybrid 
bands occur to us. The two bands may be thought 
to originate from two entirely different vibra- 
tional transitions. As they both have to be anti- 
symmetric OH stretching frequencies, one might 
be the third harmonic of the antisymmetric 
frequency, and the other the combination fre- 
quency of the antisymmetric stretching frequency 
and the second harmonic of the symmetric 
stretching frequency. This explanation seems 
improbable, however, since the expected fre- 
quency difference is much larger than the 
observed doublet separation. 

The two bands may be imagined to be members 
of a sequence of bands all of which involve the 
same vibrational transition but which involve 
absorption from different excited states. Thus 
the strongest band may originate in a transition 
from the ground state while the band displaced 
7.40 cm to lower frequencies originates in a 
vibrational transition that is different in that the 
torsion oscillation is excited in the lower and 
upper levels. In this case, however, there should 
appear a third member of the band sequence 
which is displaced about 7.40 cm from the low 
frequency band and which is about two-thirds 
as strong as this band. Actually one is unable to 
discern a trace of a third member of the imagined 
band sequence; consequently we have decided 
that the two observed bands are not members of 
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a sequence of bands involving absorption from 
different excited states of the torsion vibration. 

According to theory the vibrational levels of 
non-planar hydrogen peroxide are doubled. There 
will be a sizable doublet splitting if one of the 
barriers of the potential hindering the internal 
rotation of the OH groups is low. We are inclined 
to believe that this situation most plausibly 
accounts for the two bands of nearly equal 
intensity and identical rotational structure; 
hence we have tentatively decided that the 
bands arise from transitions between well-split 
doublet levels. The separation of 7.40 cm of the 
observed bands probably represents the difference 
in the splitting of the ground state and the upper 
vibrational state. The doublet splitting which 
results from the double minimum and the one 
low barrier in the torsion potential may be 
expected to be different in the ground state from 
that in the excited vibrational state since there 
should be a rather strong interaction between the 
OH stretching motion and the torsional motion. 
Stretching of the OH bonds changes the moments 
of inertia of the torsionally oscillating OH 
groups and also changes the hindering potential. 

At present we are unable to give an explanation 
of the apparent quadruplet character of part of 
the observed rotational structure. 

In order that our decision concerning the 
origin of the bands may be confirmed a further 
spectroscopic study of hydrogen peroxide should 
be made. We hope to undertake soon the study 
of the fundamental bands and thereby obtain 
some additional information regarding the spec- 
troscopic behavior and the nature of hydrogen 
peroxide. 

In conclusion the authors would like to express 
their appreciation to Professor Richard M. 
Badger for his helpful advice and continued 
interest in this study. 
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The simple London theory of the diamagnetic anisotropy of organic molecules containing ben- 
zene rings is modified to include the effect of a non-orthogonality of 25 percent between atomic 
wave functions on adjacent atoms. Except in the case of diphenyl the effect of such a modifi- 
cation is surprisingly small ; that is, although the contributions to the anisotropy from individual 
electronic orbits is materially changed the resultant contribution of all the electrons is hardly 
altered. The significance of this result is discussed. 





INTRODUCTION 


ECENTLY a theory has been developed by 
F. London! to explain the large diamagnetic 
anisotropy of certain aromatic molecules. Con- 
sidering the crudeness of the model the results 
are in remarkably good agreement with experi- 
ment, especially in view of the fact that the 
diamagnetism is far more sensitive than the 
energy to the form of the approximate wave 
functions. The purpose of the present note is to 
examine the effects of a certain improvement in 
the model, secured by taking into consideration 
the large non-orthogonality between the atomic 
wave functions. 


GENERAL FORMULAE 


The present discussion follows closely the 
notation and procedure of London’s paper.! The 
method used is the ‘‘approximation of tight 
binding’ as applied to aromatic molecules, for 
example, in the work of E. Hueckel.? Each 
atomic function satisfies the wave equation in a 
magnetic field with the origin of the vector 
potential chosen at the given atom, but in 
setting up the molecular orbitals the atomic 
functions are subject to a gauge transformation 
which transfers the vector potential to which 
they refer to a common origin at the center of the 
molecule. To make the discussion definite we 
consider the secular equation for naphthalene.* 
When this equation is revised to take account of 
the overlap integral 6 between adjacent atomic 


1F. London, J. de phys. et rad. 8, 397 (1937), designated 
as I. H. Brooks, J. Chem. Phys. 8, 939 (1940). 

* E. Hueckel, Zeits. f. Physik 70, 204 (1931). 

* Reference 1, p. 403. 


wave functions, it has the form: 
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where 


m=1—6(x+a), x=(E-—Eo)/Wi, a=E,/W,, 


Ey= f ux Vurdr, Wi f ux Vuxssdr, 


fi2=(e/he)3H (x1y2—X241). 


(xx, ¥x coordinates of atom K). 


€:=exp (27if12), 


If we divide each row or column of this determi- 
nant by 1—6(x+a), we obtain a secular equation 
which is formally identical with London’s, except 
that the energy variable is replaced by 


x 


C9 renee, 3 
. 1—6(x+a) (3) 


This secular equation we shall hereafter designate 
by P. Let y; be a root of P; then from (3) we 
obtain immediately that 


x= (1—da)(1+591)-'91 (4) 
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TABLE I. Comparison of anisotropy with and without non- 
orthogonality. 








y” ” fad 2” 


Benzene 
0 —0.0278 —0.0123 
1/(24/3) —0.0139 —0.0195 
—1/(2/3) —0.01389 —0.0195 
Sum=—0.0555 —0.0513 


Naphthalene 
—0.0322 
0.1382 
—0.1985 
—0.3615 
0.3333 


Sum =—0.1207 








—0.033 
—0.177 











Sum =—0.192 


Sum = —0.123 











Thus the eigenvalues of Eq. (1) are immediately 
expressible in terms of the eigenvalues of P, 
which have already been obtained by London in 
I, provided we make the usual approximation of 
disregarding integrals between all but nearest 
neighbors. From Eq. (4’) we see that (W;—6Eo) 
actually plays the same role as W, in London’s 
theory, i.e., it is the analog of the ‘‘hopping”’ 
integral. The denominator 1+éy; is just a 
renormalization factor. 

In the presence of a magnetic field the 
eigenvalues y; depend in general on the flux 
2nf = 2xeS/hc, and so also does x, so that we have 


y=y' ty" (2af)+y"" (2xf)?*, 
x=x' +x" (2rf) +x" (2rf)?. 


Making these substitutions in Eq. (4) and 
equating the coefficients of successive powers of 
(2xf), we obtain the following results 


x’ =(1—6a)(1+édy’)-y’, 
x’ = (1—da)(1+dy’)-*y”, 
1—éa 1—éa 
mr tt 9 
wai , - a , 7 
(1+6y’)? (1+6y’)8 


For the anisotropic part of the susceptibility it 
is readily shown that 


Xp= —2| Wi| x!” (2xeS/hc)?, 


(5) 


(6a) 
(6b) 





(6c) 


(7) 


HARVEY BROOKS 


where x, is the contribution of a single electron 
and the summation is carried out over all the 
occupied orbitals, x’”’ being taken from Eq. (6c). 

This general theory is applicable to each of the 
six aromatic compounds considered by London. 
In formula (6c) the quantities y’, y’’, and y’”’ have 
already been obtained in London’s paper, since 
they are simply the roots of a secular equation of 
type P, which is identical with London’s. 
Formula (6c) is a good illustration of the state- 
ment in the introduction that changes in the 
wave function affect the diamagnetism much 
more significantly than the energy, for the 
renormalization factor occurs squared and cubed 
in the field dependent parts of the energy. 


NUMERICAL RESULTS 


In Table I the calculations have been made 
with formula (7), assuming 6=}.‘ All the nu- 
merical values given are measured in units of 
(1—é8a). The column headed y’ gives the non- 
magnetic part of the energy of the given orbital, 
while the columns headed y’” and x’” are 
proportional to the contributions to the ani- 
sotropy from this orbital without and with non- 
orthogonality taken into account, respectively. 
Table II shows the anisotropies of the various 
aromatic molecules relative to benzene as de- 
termined with and without non-orthogonality. 

The notable feature of these results is the 
small extent to which even quite a large non- 
orthogonality alters the total anisotropy even 
though the contributions from the individual 
electrons are materially changed. The alteration 
of the relative anisotropies of different molecules 
is even smaller. The satisfactory result is prob- 
ably to be ascribed in all cases to the fact that 


TABLE II. Relative anisotropies with and without non- 
orthogonality. 








(Ax/ Axpenz) 
LONDON CORRECTED 


2.19 2.23 
3.45 3.54 
3.19 3.06 
4.46 4.75 
2.21 3.42 


EXPERIMENT 


2.11 
3.38 
3.07 
4.31 
2.20 





Naphthalene 
Anthracene 
Phenanthrene 
Pyrene 
Diphenyl 








4A. L. Sklar and M. Goeppert Mayer, J. Chem. Phys. 
6, 645 (1938). 
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we are dealing in effect with a half-filled shell, for 
which the properties of the whole half-shell may 
often be expected to be invariant of the particular 
form of the individual wave functions. The one 
exception to this rule seems to be diphenyl, for 
which the simple and corrected theories give 
widely divergent results, experiment being defi- 
nitely in favor of the simpler theory. 

Such refinements of the simple molecular 
orbital treatment should not be taken too seri- 
ously, but they do serve to illustrate the very 
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provisional and crude nature of the approxi- 
mations used in describing the structure of 
polyatomic molecules. The neglected effects are 
by no means small. The surprising thing is, 
indeed, that the crude theories give such good 
results. In the case of magnetic anisotropy the 
good results of the London theory, as pointed out 
by London, are probably to be ascribed to geo- 
metrical considerations which are invariant of 
the model, but which are well represented in the 
crude molecular orbital picture. 
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Tables of Thermodynamic Functions of Paramagnetic Substances 
and Harmonic Oscillators 


J. R. anp R. A. HULL 
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Numerical values of the partition function, magnetic moment, entropy, and temperature 
derivative of the entropy are tabulated for free paramagnetic ions with j= 4, 34, and 54. These 
tables are needed in connection with experiments on adiabatic demagnetization. As a by- 
product, the thermodynamic functions for an Einstein oscillator are computed numerically to a 
higher degree of accuracy than in the tables available in the literature. 


INTRODUCTION 


HE calculations shown in this paper have as 
their principal object the evaluation of the 
entropy and its derivative for ideal paramagnetic 
substances in terms of H/T. Such theoretical 
values are necessary in interpreting the results of 
magnetic experiments, particularly at low tem- 
peratures, and these applications are dealt with 
more fully below. 

The starting point for the calculation of 
entropy S is the partition function Q, in this case 
Le. The summation over m is from —j to +j 
where j is the angular momentum quantum 
number appropriate to the substance considered, 
and x is proportional to H/T. S and its derivative 
can then be obtained in terms of Q and its 
derivative, as is done in the first half of this 
paper. It is also possible to obtain S in a more 
indirect (though actually briefer) way; for Q can 


be rewritten as the ratio of the partition functions 
for two Einstein oscillators, and so S and its 
derivative are obtained as the differences be- 
tween the corresponding quantities for these 
Einstein oscillators. This method is used in the 
second part of the paper. 

It was felt advisable to do the calculations by 
both methods as a check on accuracy. Inci- 
dentally several other useful thermodynamic 
quantities were evaluated as intermediate steps, 
and these have also been tabulated. They include 
either explicitly or else implicitly (as quantities 
obtainable simply by addition) the magnetiza- 
tion, total energy and free energy of a paramag- 
netic substance, and the total energy and free 
energy of an Einstein oscillator. Some of these 
quantities have already been tabulated in the 
literature, e.g., the Einstein functions in Landolt- 
Bornstein, but not in sufficient accuracy to be 
always useful. 
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TABLE I. 








M/Nu (M/Nu)? = S/R (—2dS/dz)/R 


(M/Nu)? = S/R (—2d8/dz)/R 





j=1/2 
0 0.69315 
0.01000 
0.02000 


ze52e35 


SCHOAFrNwOaD 


SANS 

PNNNNNNN 
SSkerss 
RIRERSH= 


0.43923 
0.43506 
0.42366 
0.40660 
0.38535 
0.36128 
0.33553 
0.31304 
0.28260 
0.25676 
0.18653 


0.13054 


Comm SDR SSONS tia tete SS OND iim tm OS ONG Tih SSRIS ie Ww br nreRoSeS 


0 
0. 
0: 
0: 
0: 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 
0. 

0. 
0. 

0. 

0. 
0. 

0. 
0. 

0. 

1 

1 

1 

1 

1, 
1 

1 

1 

1 

1, 
2 

2 

2. 
2. 
2. 
2. 
2. 
2. 
2. 
2. 
3. 
3. 
3. 
3. 
3. 
3. 
3. 
3. 
3. 
3. 
4, 
4. 
4. 
4. 
4. 
5: 
5. 
5. 
5. 
5. 
6. 


.02053 
0.01731 0.08879 





1.38629 


0 

0.00050 
0.00200 
0: 


0.74315 
0.73874 
0.72993 
0.71763 
0.70257 


0.62576 


0.49539 


SOM DS DWM SSONG mW SS DWUSDim WIE OSHNSUR WE SSRISUBRYSwMePWRE EGER SSSes 


D> 88 8 1 Gt ee 1 ee me 69 69 09 G9 G9 9 SO OI CO CODON DDD NNN NE Sm mm mmm OSS SSS SSSsssessss2s2e22“2° 


0.01739 





I. Direct PARTITION FUNCTION METHOD 


The partition function Q has already been 
defined as 


j 


Zam, 


m=—j 


where x=guH/kT, 


g is the Landé splitting factor, and y is the Bohr 
magneton. The thermodynamic functions for a 
gram ion are obtained from Q as follows: the free 
energy F= —RT InQ; the entropy S= —d0F/dT, 
which is identical with R(InQ—Q-'xdQ/dx); 
the magnetization M=RT0a(InQ)/dH, which is 


identical with Nu2Q-'dQ/dx; the derivative 
(—xdS/dx). It should be noted that (—xdS/dx) 
= T7T(0S/8T)”4=H(dS/dH)r, where T(dS/dT)x 
is, of course, the specific heat in a magnetic field; 
and it is sometimes convenient to use the fact 
that d(S/R)/dx=43xd(M/Nu)/dx. The _ total 
energy U can be obtained immediately from the 
relation U= F+TS. The tables of these functions 
have been arranged separately according to the 
value of j (see Table I). As regards numerical 
values, g has the value 2 for the paramagnetic 
substances most commonly used, which depend 
on the paramagnetism of the free electron spin in 
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Ti, Cr, Mn and Fe ions. The values of the other 
necessary constants according to the most reli- 
able data at present available are: 


N=6.024 
k=1.380 
R=8.302 
pw=9.28 


1028 

10-'* erg/deg. 
10’ erg/deg. 
10-*! erg/gauss. 


From these values it follows that if H is expressed 
in kilogauss, then x =0.134;H/T. 


Uses 


It may be of interest to mention some of the 
uses to which these tables can be put, particu- 
larly in connection with the demagnetization 
technique. The value of M, for instance, is used 
when computing the reduction in the effective 
applied magnetic field caused by the demag- 
netizing effect of the substance in question. The 
value of S in terms of H for a given T is required 
for the interpretation of the resulting low temper- 
ature obtained by adiabatic demagnetization 
from various values of H; for the entropy at this 
low temperature in field zero is the same as the 
value of S in field H at the initial temperature T. 
The quantity (—xdS/dx), being in the nature of 
a specific heat, can be used to interpret the 
experimentally determined specific heat. 

It is also possible to use the tables to get some 
idea of the variation with temperature of the 
entropy and specific heat in zero field of an 
actual paramagnetic substance, at temperatures 
so low that it is no longer ideal. Namely, one may 
assume as a rough approximation, as did Kurti 
and Simon,! that the ground state of such a 
substance is split into 2j7+1 equidistant levels 
separated by an energy k@, where @ is a charac- 
teristic temperature for the substance. The 
entropy and specific heat are thus given by 
putting x =6/T in the tables of Sand (—xdS/dx). 
By using the approximation formulae given 
later, it is seen that at temperatures high com- 
pared with @ (i.e., x small), the specific heat is 
inversely proportional to 7? and the entropy 
defect, reckoned from the entropy when x is 
zero, is also inversely proportional to 7°. 


1N. Kurti and F. Simon, Proc. Roy. Soc. 149, 152 
(1935). 


Range of values 


These tables were calculated for their utili- 
tarian value in work with the magnetic method. 
The values of j were therefore confined to 14 (for 
Tit++), 3¢ (for Cr+++), and 5¢ (for Mn++ and 
Fett++). This choice includes the substances 
most frequently used with this technique. At 
some later date'* it may be possible to extend 
them to 74 for Gd (whose salts are not likely to 


TABLE I.—Continued. 








InQ M/Nu 
j=5/2 


(M/Nu)? 8/R (—2dS/dz)/R 





1.79176 1.79176 
1.79118 
1.78943 
1.78653 
1.78248 
1.77731 
1.77103 


0 
0.11664 


San oBeOnwmnornwonanm 


0.83099 
0.84557 
0.84941 
0.84536 
0.83563 
0.82184 
0.80519 
0.78648 
0.76630 
0.74510 
0.72318 


SRR SODO> RWS DONG WIR SS DONS RWS OS DNS Rw SSIS ine wewiinbm ere SSeS 


SP SY OS C7 oe om eo G0 G0 G9 GO GD GO DO DO BW DO I ON NN NNN ND Se ee ee eee OOD ODS OS SOO SSS SSSesSssSssSso 


988800 
3277100 








18 Note added in proof.—A table of the magnetization 
and entropy appropriate to j7=7/2 has just been published 
by Giauque in the February, 1941 issue of the Journal of 
the American Chemical Society. 
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TABLE II. Einstein functions. 
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In(e* —1) Sz/R (e# —1) 
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0.16886 

0.15719 A . 0.49627 

0.14624 / . 
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0.12638 . 

0.11739 . . 0.41290 

0.10898 : . 

0.10113 ° e 0.37429 

0.09380 A , 

0.08695 . 5 0.33799 

0.08058 : . 

0.07463 5 J 0.30408 

0.06394 : d 0.27264 

0.05469 : ld 0.24364 

0.04671 - r 

0.03983 . x 0.19277 

0.03392 5 J 

0.02885 . J 0.15083 

0.02450 ° A 0.13292 

0.02079 : y 0.11683 

0.01761 . \ 

0.01491 ‘ - 0.08967 

0.01261 r 

0.01065 ; 0.06829 

0.00899 . 0.05942 

0.00758 ; , 0.05162 

0.00639 A r 

0.00538 . ; 0.03876 

0.00453 7. J 

0.00381 ‘ J 0.02893 

0.00320 ° X 0.02495 

0.00268 A A 0.02148 
0.01587 


0.01000 
0.00624 
0.00387 
0.00238 
0.00145 
0.00088 
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be used so much in the future) and to % corre- 
sponding to certain nuclear spins.? 

Similarly the range of values of x is that likely 
to be useful in practice. For M and S, the range of 
x extends up to 6.0, corresponding to fields of 44 
kilogauss at 1°K; for still higher values of x the 
approximation formulae given below will usually 
be sufficient. In the tables of (—xdS/dx), x 
extends only to 2.0, as in this work the calcu- 
lations involving (—xdS/dx) are usually confined 
to small fields. For 7=44 and 34, however, the 
tables show (—xdS/dx) at some values of x 
higher than 2. These have been obtained from 
the table of Einstein functions (see Table IT). 
For example, (—2dSz/dz) with z= 2.2 was one of 

2 The unit of nuclear magnetic moment is 1/1840 times 


the Bohr magneton; it is not possible to give any general 
rule for the nuclear g value. 


the quantities calculated to check (—xdS/dx) for 
j= with x=1.1; it is now used as one of the 
quantities to give (—xdS/dx) for j= with 
x=2.2.33 
APPROXIMATION FORMULAE 
For small x one has the developments 
M/Nug = 3xj(j+1) 
—(1/720)[(2j7+1)4—1]x?+---, (1) 
S/R=1n(2j+1) —§x°7(G+1) 
+(1/960)[(2j+1)4—1]}x4+---. (2) 


3 If still further values of (—xdS/dx) are required, to an 
accuracy greater than that obtainable by graphical differ- 
entiation of S, the Einstein function method given later 1s 
by far the shortest and most convenient. 





Qr= 


The 
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For large x one may use the relations 
M/Nug=j-—e*+e-*, (3) 
S/R=(x+1)e-*F (x+3)e7". (4) 


Here the upper or lower sign choice is to be 
employed, according as j= or j21. 

If an accuracy of 1 percent is desired it suffices 
to use merely the first term of (1) or (4), and two 
terms of (2) or (3), provided x has the following 
range of values: <0.32, 0.24, 0.12 in (1), >2.8, 
2.2, 2.0 in (3), <0.80, 0.40, 0.32 in (2) and 
>4.8, 4.4, 4.4 in (4). Here the first, second, and 
third choices apply according as 7 = 4, 34, 54. The 
first term of the formulas for xdS/dx obtained by 
differentiating (2) and (4) is enough to give one- 
percent accuracy if x<0.20, 0.12, 0.06 or 
x>6.0, 5.4, 5.4. 


II. EINSTEIN FUNCTION METHOD 


For an Einstein oscillator of frequency v the 
thermodynamic functions are again derived from 
the appropriate partition function; the free 
energy and total energy are, however, slightly 
different according as the energy of the lowest 
quantum state is taken as being hv/2 or zero. In 
the following treatment we consider always N 
oscillators. If the energy states are hy, 3ghv--- 
etc., then the partition function is 


Op =e /2kT 4 g—Bhv/ RT 4 g—Shv/ kT... 
=e/2/(e*-—1) where z=hv/kT. 
The free energy is 
F=RT-[In(e*—1) —32], 


while the total energy is 


4 z 
u=Rr|-+ | 
2 e-—1 


If the energy states are 0, hv, 2hv, --- etc., i.e., if 
there is no zero-point energy, then the expres- 
sions for the free and total energy are obviously 
diminished by %hv as compared with the 
above values. The entropy Sg is given by 
R(InQ—Q-xdQ/dx], and is thus R[z+2(e?—1)— 
—In(e?—1)] regardless of whether there is a zero- 
point energy and (—z2dSz/dz) is Re’e?/(e*—1)?. 
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The accompanying tables show Sg/R and 
(—2dSz/dz)/R, also z/(e?—1) and In(e*—1) so 
that the free energy and total energy are easily 
obtained for either partition function. 

The importance of the Einstein functions for 
our particular purpose lies in the fact that the 
partition function 


? 
bens 
m=—j 

for a paramagnetic substance is the ratio of two 
Einstein functions e*/?/(e7—1), where z has the 
values x and (2j+1)x, respectively. Since S is 
given by R[InQ—xd(InQ) /dx], it isseen that S/R 
for the paramagnetic substance is obtained as the 
difference between the two quantities Sg/R for 
Einstein oscillators with z equal to x and (2j7+1)x, 
respectively. In other words 


S=Sz(2=x)—Se (2=2jx+x). 


Analogous relations apply to —xdS/dx. This 
useful property has been pointed out by Kurti 
and Simon. 

The selection of values of z in the tables, which 
seems somewhat arbitrary at the higher values, is 
just that required to cover the same range of S 
and (—xdS/dx) as in the direct partition function 
method. The resulting values of Sand (—xdS/dx) 
have not been tabulated again, as they are 
identical with those already obtained. As already 
mentioned, this agreement serves as a par- 
tial check on the numerical accuracy of our 
calculations. 


ACCURACY 


While making the calculations six significant 
figures were used, and the results have since been 
rounded off to five figures. The calculations were 
made as far as possible on an electrical machine, 
which accounts for this large number of signifi- 
cant figures compared with the intervals in the 
series of x and z values. There may, however, be 
occasions when it will be useful to have available 
the five significant figures, and it is hoped that 
the error is nowhere greater than one in the last 
figure. 

We wish to record our indebtedness to Dr. N. 
Kurti for many helpful discussions. 
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The ordinary formula for the rotational energy of a simple rotating vibrator, when applied to 
a diatomic molecule, consists of the first few terms of a series that becomes poorly convergent for 
levels of high rotational energy. Inversion of the series results in a better approximation. An 
illustration is given using the results of calculations for the normal state of hydrogen. The 
inverted series is especially suited to derivation of approximation formulas for the partition 
function. Such an approximation formula should not be carried to powers of T so high that lack 
of coefficients for the ordinary formula introduces significant error. 





N the course of some calculations of thermo- 
dynamic properties of hydrogen! from spectro- 
scopic data, a difficulty was encountered arising 
from the form of the expression for the rotational 
energy of a diatomic molecule. In ordinary 
spectroscopic work with heavy diatomic mole- 
cules the formula 


is fully adequate since the second term is very 
small. For hydrogen and for large energies this 
is no longer the case. A third term is in use in 
spectroscopic work on hydrogen and, to extend 
the application to higher energies, a fourth can 
be evaluated using relations from the theory of 
the simple rotating vibrator. The resulting 
formula for the rotational energy, 


E,=B,K(K+1)+D,K2(K+1)? 
+F,K#(K+1)!+H,K(K+1)4, (2) 


while adequate for a somewhat greater range, 
is still not satisfactory since it is the beginning 
of a series that is too slowly convergent or even 
divergent for large but possible values of rota- 
tional quantum number, K. The formula is, in 
fact, not very suitable for calculating thermo- 
dynamic properties of hydrogen above 2000 
degrees K. Some other treatment is therefore 
necessary. 

Using a modification of the Klein-Rydberg 
method,? the molecular potential energy for the 
ground state of hydrogen was determined as a 


1 A project of the National Research Council Committee 
on Thermodynamic Data for Chemical Industries, Dr. 
F. R. Bichowsky, Chairman. 

2 R. Rydberg, Zeits. f. Physik 73, 376 (1931); O. Klein, 
Zeits. f. Physik 76, 226 (1932). 


function of internuclear distance and with this 
the quantum integral was evaluated numerically 
for several combinations of large rotational 
quantum number and total energy. The results 
indicated that, at least for the particular spectro- 
scopic constants employed, the formula became 
adequate for the purpose when continued indefi- 
nitely as a geometric series with the ratio given 
by the last two terms of formula (2) and with 
the added part summed with the usual formula. 
Figure 1 shows, for v=0, the estimates which 
would be made for E, as a function of K, using: 
for curve 1, the first term of formula (2); for 
curve 2, the first two terms; for curve 3, the 
first three terms; for curve 4, the first four terms; 
for curve 5, the sum including the continuation 
as a geometric series; and for curve 6, the 
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Fic. 1. Various estimates for the rotational energy of 
the hydrogen molecule as a function of rotational quantum 
number for the lowest vibrational state. Curves 1, 2, 
and 4: Formula (2) evaluated to 1, 2, 3 and 4 terms, 
respectively. Curve 5: 4 terms of formula (2) with a 
geometric series continuation. Curve 6: Result of detailed 
numerical investigation. Curve 7: Formula (3) evaluated 
to 4 terms. 
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interpolated results of the numerical evaluation 
of the quantum integral previously mentioned. 
Several other usable methods of representing 
the relation between E, and K have been 
suggested to the writer. One of the most inter- 
esting of these is the inversion of the series of 
formula (2). This leads to 
1 D, (2D,?—B,F,) 
K(K+1)=—E,—-—E,?+ E,? 
B, B,? B,§ 
Ph ssninsiiare iinet 


B,’ 








Eyt+--+. (3) 


In using this, the total rotational and vibrational 
energy is to be G,+£, where G, is the purely 
vibrational energy. In contrast to formula (2), 
in which the terms alternate in sign, formula (3) 
as far as evaluated is rapidly convergent, with 
all terms positive, for all values which £, can 
properly assume. In Fig. 1, curve 7 shows the 
result given by the first four terms of formula (3) 
for v=0. Comparison of this with curve 4 shows 
that, without including any estimates of re- 





Kmin 


= {exp ( 


Continued integration by parts gives for this 


- [exp (=) +(-) 


If (E,)max is great enough so that 


*@*[K(K+1)] 
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mainders, formula (3) is definitely superior to 
formula (2) for large values of K and E,. 

Formula (3) is particularly adapted to use in 
deriving one type of approximation formula for 
the partition function.’ In the customary deriva- 
tion, the summation 


DK {exp (—c2(G,+£,)/T)}(2K+1) 


is represented in accordance with the Euler- 
Maclaurin summation formula as the sum of the 
integral 


Kmax 
f exp (—¢2(Ge+E,)/T)} (2K +1)dK 
K 


min 


and other terms, known in this application as 
the Mulholland correction,‘ and obtained from 
the integrand and its derivatives. In treating 
the integral it is preferable not to break up 
part of 


exp (—c2E,/T) 


into a power series expansion in K(K+1) but 
rather to use the transformations: 


nr xp (- et) | oK+ax= fo fo (SE) fat +0 


‘=Kmin 





2Uy E v)max Ey d| K(K+1 
=)I f' ) lex (— )| [K(K+ dn. 
T T dE, 


(Ey)min 








aS 


dE,? Cs dE, 
i 4d‘(K(K+1) ] (EZ ») max 
2e- 


dE,4 (EB »)min 





C2 


exp (—c2(G,+£,)/T) 


makes the contribution of the upper limit negligible and if (E,)min=0, this becomes, by formula (3), 


6B,F, T*120B,D,F,— Rntenentailnn 








—) + T*12D3— 


oo YEE 


oa B,' 


(4) 


* W. F. Giauque and R. Overstreet, J. Am. Chem. Soc. 54, 1731 (1932); L. S. Kassel, J. Chem. Phys. 1, 576 (1933); 
H. L. Johnston and C. O. Davis, J. Am. Chem. Soc. 56, 271 (1934). 


* H. P. Mulholland, Proc. Camb. Phil. Soc. 24, 280 (1928). 
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It must be pointed out that evaluation of formula (4) cannot be extended to higher powers of T 
without first evaluating additional coefficients of formula (2).5 For example, the coefficient of T‘ 
in formula (4) is thrown off by several hundred percent by setting H, =0. In addition to the numerical 
error which would be introduced, this could lead to the inference that the series is somewhat less 
convergent and hence less adapted for calculation than is in fact the case. 

With g representing the average statistical weight due to nuclear spin degeneracy and with the 
zero of the energy scale chosen so that Gp=0, the internal partition function is the sum over v of 


co(G,+E,) 2G, 7.4 2D, T3 
a oo ae | al ina er) yar er are 


12D,?-—6B,F, T7*120B,D,F,—24B,?H,—120D,' 1 1 cB, 
x — iceadiinel 
B,§ Cot B,' 








4 C2?B,? 


$18 fF HS F* 

A few terms of the Mulholland correction have been included in the formula. They are to be 
used with caution since the true correction is not properly developable in powers of J-! and the 
apparent convergence is too optimistic an indication of the closeness of the approximation, especially 
with T small. As is well known, this correction does not distinguish between sums over even and 
over odd values of K. 

It is desired to acknowledge the encouragement by Dr. F. G. Brickwedde of the calculation of 
thermodynamic properties of hydrogen at high temperatures from spectroscopic data. Several 
suggestions made by Professor David M. Dennison have been very helpful in leading to this paper. 


5 For the simple rotating vibrator, see J. L. Dunham, Phys. Rev. 41, 721 (1932); I. Sandeman, Proc. Roy. Soc. Edin- 
burgh 60, 210 (1940). 
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The condition which has to be fulfilled in order that the Bose-Einstein condensation be 
accompanied by a discontinuity of the heat capacity at the condensation temperature is dis- 
cussed. The nature of the transition depends on the distribution of the energy levels of the 
individual particles forming the gas. The condensation is studied for several types of individual 
energy spectra. It is shown that all those Bose-Einstein (B.E.) assemblies which admit a non- 
ambiguous definition of the statistical distribution function present, beside interesting analogies, 
a qualitatively different behavior from that observed in liquid helium. It appears that as long 
as the transition in momentum space of non-ideal B.E. systems is not understood, the possibility 
of connection between B.E. condensation and the transition observed in liquid helium remains 
an open problem. 


N a previous paper! we have discussed the 

condensation of B.E. ideal gases in presence 

of conservative force fields which acted directly 

on the motion of the center of gravity of the gas 

atoms. It has been assumed there that the 

classical Hamilton function could be used for the 
1L. Goldstein, J. Chem. Phys. 9, 273 (1941). 


description of the motion of the individual gas 
atoms. This was equivalent to the assumption 
that practically no serious error was introduced 
by leaving aside the quantization of the motion 
of the individual atoms in those fields. The 
modification brought about by the field in the 
properties of the gas around the condensation 
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temperature concerned uniquely this last quan- 
tity. The presence of the field increases the 
transition temperature above that of the free gas 
with the same concentration. The condensation 
remained a phase transition of the third order, 
since energy and heat capacity at constant 
volume were continuous at the transition temper- 
ature, the discontinuity presenting itself was that 
of the derivative, with respect to the tempera- 
ture, of the heat capacity. Apart from the 
intrinsic interest in the study of the peculiar B.E. 
condensation under different conditions, it 
seemed that the eventual possibilities of con- 
nection between this effect and the well-known 
transition effect occurring in liquid helium? 
justify a more general discussion of some of the 
properties in thermal equilibrium of B.E. gases 
than that made on a classical basis. We will be 
concerned here mainly with the conditions which 
assure the reduction of the order of B.E. transi- 
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tions and introduce, therefore, a discontinuity in 
the heat capacity at constant volume of the 
system at the transition temperature. Such B.E. 
systems presenting a transition of the second 
order seem to provide quite an attractive picture 
of what may eventually happen in liquid helium. 
However, beside the apparent analogies between 
B.E. condensation and some features of the 
transition in liquid helium there is also an im- 
portant difference of a qualitative nature which 
will be discussed below. Evidently, no secure 
model of liquid helium can be based on ideal B.E. 
systems or B.E. systems under appropriate fields, 
self-consistent fields included. The internal field 
of non-ideal B.E. assemblies might modify con- 
siderably the transition in phase space. The study 
of the eventual connection between the B.E. 
transition and that observed in liquid helium 
should be based on the analysis of the thermal 
properties of non-ideal B.E. fluids. 


We have to recall briefly some properties of the ideal B.E. gas in absence of any force field. 
The average energies, at temperatures 7, below and above the transition temperature 79, are given, 


respectively, k denoting Boltzmann’s constant, by 


E_=3NkT(T/T»)'G(0)/F(0), 


E,=4NkTG(a)/Fla); F(x) =¥ e/a; 


(1) 
G(x) = e#/A8!2; (2) 


A=1 


where the always positive fundamental parameter a@ is defined by the total number N of spinless 
atoms of mass m in the volume V through the equation, 


N=(V/h*)(2amkT)'F(a). (3) 


The heat capacities at constant volume are, respectively, 


15 
' LEE NE Ee (4) 


3 d da 
dE... /dT =—Nk\ G(a)/ F(a) +7T7—(G/F)—l. 
aT = [G1a/r(a) +7 (G/F) | 


Since, from (3), 


(5) 


da/dT'=—>-Fla)/F(a); Pla)= = —5 e/a! (6) 
a/ OT (a)/F'(a); a at a ¢ 


one gets 
Ey 


dT 


a 


d 3. 75 3 
—*=-wil -G(a)/Fla)+-Fla)/ F(a) | (7) 


At high temperatures the series G(a), F(a) and F’(a) reduce, practically, to their first term e~* and 
* F. London, Nature 141, 643 (1938); Phys. Rev. 54, 947 (1938). 
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—e~*, respectively, because a is large. The heat capacity decreases and tends asymptotically toward 
the classical value 1.5Nk. At low temperatures, just above the transition temperature, a~0, and the 
first term on the right-hand side of (7) tends toward 15/4NkG(0)/F(0) as does (dE_/dT)ro. The 
second term on the right-hand side of (7) tends toward zero since the series F’(0) or 


—-> 1/M 
A=1 








diverges. Therefore, following (6), 


(da/dT)T)=0; 







(da/dT)r>T.2 0. 








It is instructive to write dE,/dT in the following form 


dE, 15 3. da 
—— =—NkG(a)/ F(a) -—-NkT—, (9) 
dT 4 2 dT 





which shows clearly, with (4), that the continuity of the heat capacity at constant volume of ideal 
B.E. gases at the transition temperature is due to the fact that the derivative of the fundamental 
parameter a with respect to the temperature vanishes at the transition point. 

The second derivative of the average energy of the gas is, below 7», 


@E_ 45 Nk/T — - 
ere oe, 60/0), 














which is always positive. The heat capacity at constant volume increases until, approximately, 
1.92Nk at Ty. Above the transition temperature the heat capacity decreases monotonously from the 
preceding value and tends to 1.5Nk. An elementary calculation gives 

@E, 45 NkG(a) 9 Nk Fla) 27 Nk[F(a)]}?F’(a) d°F 


ips ni aie tiieds aii: ipemceslions ta ; F’=—. (11) 
dT? 8 T F(a) 4T F(a) 8 T [ F’(a) }? da? 


























At high temperatures G(a), F(a) and F’’(a) reduce practically to e~*, F’(a) to —e~* and lim? targe 
@E,/dT? vanishes. The discontinuity of the derivative, with respect to the temperature, of the heat 
capacity at constant volume at the transition temperature is, with (10) and (11), 


@E @E, d*E_ 
dT*/ 17 \dT? / 7 dT? / 1% 
27. Nk (F(a) F(a) 
=— — lim 
8 T 9° [F'(a) }? 
















(12) 








which is negative. 

















Our main problem is to find under which con- parameter is finite at that temperature. It is 
ditions the phase transition of the third order as therefore necessary to investigate the B.E. gas 
defined by (12) is reduced toa phase transition of _ under such conditions which assure a finite value 
the second order, i.e., with a discontinuity of the for da/dT at the transition temperature. 









heat capacity at constant volume at the transi- We must remember here that a is defined at 

tion temperature. temperature T by the total number N of the gas 
It is clear from what precedes that the B.E. atoms through 

ideal gas cannot present a discontinuity of the N=T, 1%; 

heat capacity at its condensation temperature as (13) 






unless the derivative da/dT of the fundamental = Di gi/exp (a+ 8W;)—1, 
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where n; is the average number of atoms in the 
state of energy Wi, g; is the weight of this state, 
Bis 1/kT and the summation extends over all the 
states W; of the energy spectrum of the indi- 
vidual atoms. It will be assumed throughout this 
paper that under all circumstances the possibility 
of the non-ambiguous definition of the 1; or the 
distribution function ;/g; is maintained. There- 
fore, the only way to modify the function 
a(N,T) is to introduce or impose new energy 
spectra (W;, gi) for the motion of the individual 
gas atoms. This can be achieved by assuming the 
application of external fields or by allowing the 
existence of self-consistent fields arising from the 
interaction of the atoms of the system. In any of 
these imposed fields, the motion of the individual 
atoms and therefore their energy spectrum will 
be supposed to be known completely. We will 
discuss now several types of individual energy 


spectra which in some cases will assure a B.E. 


condensation accompanied by a discontinuity of 
the heat capacity at constant volume at the 
transition point. 


II. 


Consider the case where the individual atoms 
of the system have a series of / discrete stable 
levels of energy —Eo, — Fi, ---, —Ei, +++, —Ei, 
followed by a continuous spectrum of positive 
energy levels corresponding to practically free 
atoms. We are not concerned here with the 
actual existence of the type of fields which give 
rise to the preceding level spectrum or similar 
spectra and, therefore, no attempt will be made 
here to justify them. For the following discussion 
it is not indispensable that the positive energy 
levels form an entirely continuous spectrum. 
This may have a series of levels and bands which 
for sufficiently high energies would be followed 
necessarily by the free-particle continuous spec- 
trum. The inclusion of such positive energy 
levels and bands would only complicate the 
calculations without increasing the generality of 
the discussion. Another simplification can be 
made without impairing the generality of the 
results, namely, that the discrete levels will be 
supposed to be nondegenerate and their number / 
to increase eventually in such a way that 
lim 1/N—0 when N— ~. Under these conditions, 
the average number of spinless atoms of mass m 
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in one of the discrete states, say the ith, is, at 
temperature 7, with 8 denoting 1/kT, 

n:=1/exp (a—BE;)—1 (14) 


and in any one group of the continuous energy 
states with momentum between p and p+dp, 


dn(p) =g(p)dp/exp (a+Bp?/2m)—1, (15) 


g(p) being the level density in the level group in 
question. If there be N atoms in the volume V 
the parameter a is defined by 


l 1 
N= 2 
i-0 exp (a—BE,)—1 





——f p*dp 
h? Jy exp (a+6p2/2m)—1 





(16) 


Since no term of the right-hand side of (16) can 
become negative, one must have a2 BE, at any 
T. With 

a’ =a—BEo (17) 


one gets, after integration over the continuous 
spectrum and division by N, 


1/N l 1/N 
1= +2 
exp (a’)—1 i=1 exp (a’+8(E)—£,))—1 
( ~) F(a’ +BEo) 
To F(0) 








» (18) 


where 7> is the condensation temperature of the 
free gas with the same concentration N/V as the 
gas under discussion and F(x) has been defined 
by (2). At all temperatures 7 >T7 > where in the 
last term on the right-hand side a’ must be finite 
in order to reduce this term to unity or less, the 
first two terms on the right-hand side can be 
neglected on account of their extremely small 
numerators 1/N. a’ is then defined by 


BB pechennal 
7 FO) 





(19) 
This can also be written as 
V Ey 
=—(2rmeT)'#(a!+—.), (19a) 
h® kT 


which shows that at all these temperatures the 
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maximum number of atoms which the phase 
space can accommodate is 


V 
iE Pes 5. (20) 


If supplementary particles are added to the gas 
they must go to those states which are not 
included in (19a), i.e., the stable discrete levels. 
In other words, with N and V kept constant, the 
root T; of the equation 


fl x F(Eo/kT) 


=: 21 
To F(0) i 


is that temperature below which some of the 
atoms must condense or go over to those states 
which have been left out of account. It is clear 
that 7,27», a result which is similar to that 
obtained in the semi-classical treatment of the 
influence of conservative force fields on the B.E. 
condensation. Remembering the behavior of 
F(x) one sees clearly in (21) that the larger Eo/k 
the larger is T;, which is readily understandable. 
In order that condensation occur, the concen- 
tration N/V being kept constant, one must 
decrease the temperature from regions of high 
temperatures T>TZy;. At such temperatures, 
practically, no particles can be found in the 
group of stable levels: they are all empty. This 
means that kT Epo. Since Ty; is the smallest 
temperature at which the lowest levels remain 
empty, it is clear that J; has to vary in the same 
direction as Ep. 





Eo l E; 
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At temperatures T<Ty the last term on the 
right-hand side of (18) is certainly smaller than 
unity, (21) having but one root T;. The first two 
terms have to be large enough now in order to 
maintain the equality. The sum or the second 
term can still be neglected since the denominators 
of this sum are always finite whereas the numer- 
ators are extremely small, N being enormous. 
The first term becomes of the order of unity if 
a’~ N-, i.e., extremely small. It can safely be 
neglected in the argument of F(a’+Eo/kT) and 
one gets, in this approximation, 


1 (=) F(Eo/kT) 


1=—— ——_—— 
Na’ F(Eo/kT;) 


by using the definition (21) of 7;. This gives a’ 
and also the number of atoms No which at 
temperature T<Ty; are condensed in the lowest 


state: 
-1-(— -) F(E)/kT) 
T; F(Eo/kT;) 
This shows that as long as kT; is not extremely 
large in comparison to the smallest separation 6 
of the stable levels, i.e., 6/kRT;>N-, only the 
lowest level determines the condensation. The 
remaining stable states are practically always 
empty. If 6/kT; were of the order of N-, then 
instead of the one lowest state a group of 
states would play a predominant role in the 
condensation. 
One finds then easily for the average energy of 
the system 


3 V 
= +, eee a) (24) 





a » 
i=1exp (8(Eo— 
Eo U E; 


exp (a’)—1 


E))—-1> 2h 





Rin = 
+ exp (@’)—1 


AC,=C.4.—C,- 
t E;exp ((Zo— 


my 
i=1 exp (a’+8(Eo— 


3V Eo 
+-——(2xmkT)ikTG ’+—), 
E;))-—1 2h’ kT 
since no matter how small a’ is below 7; it must necessarily be taken into account in No. These 
energies are naturally continuous at T;. 
One gets, after a simple calculation, the following expression for the discontinuity in the heat 
capacity at constant volume at the transition temperature 7,, 


(25) 


— E;)/kT;) 





i=t (exp ((Zo— 


(=), -ver(F) 
E)/kT,) —1)\dT/ 1, 2 ‘\aT/1, 
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by using (21) and the expression (3) for To'F(0). That the lowest level does not contribute to AC, is 
evident since a’ and da’ /dT are continuous at 7;, as can be proved on the basis of the exact Eq. (18). 
The contribution to AC, of the stable states comes from E, where in the average energy of these 
states a’ figures explicitly. That this contribution is smaller than that of the positive energy states can 
be seen easily. Therefore AC, is different from zero unless (da’/dT)r; is itself zero. In any case 


(da’/dT) 2 0. In order to find it in the present case one has but to use (19a). One finds thus 


F(Eo/kT;) Eo 


(27) 





da’ 3 
T(—) "9 kT, 
dT Ty 2 [ —F’(Eo/kT;) ] kT; 


which is finite in any case; the first term on the right-hand side is also finite since 


—F'(x)=> e*/v; x=Eo/kT; 


A=1 


is convergent. The discontinuity AC, is finally 


S 6 


E; exp ((Eo— Ei) /kTs) 








AC,={— 2D 


This discontinuity being negative is of the same 
type as that observed in lambda-point transi- 
tions. One may say, therefore, that as soon as the 
individual level spectrum of the gas atoms includes 
at least one discrete stable level which is followed 
by the free particle continuous spectrum the B.E. 
transition becomes that of the second order. 

Beside the case just discussed, the nonvanishing 
of the derivative (da’/dt)7;, assuring thus a B.E. 
transition accompanied by a discontinuity of the 
heat capacity at the transition point, can be 
eventually achieved in a number of different 
ways corresponding to the motion of the gas 
atoms in appropriate fields. The modifications of 
the individual energy spectrum of the atoms 
which assure a B.E. transition of the second order 
are the consequences of the quantization of the 
motion of these atoms under given conditions. As 
long as the quantization is left out of account the 
B.E. transition is of the third order. However, the 
quantization does not lead necessarily to such in- 
dividual energy spectra which would assurea tran- 
sition of the second order, as it will be seen below. 

The results of the preceding discussion may 
then be summarized by saying that whereas the 
condensation phenomenon is a direct conse- 
quence of the symmetry of the system in phase 
space, the nature of the transition is determined 
by the properties of the energy spectrum of the 
individual atoms. 


F(Eo/kT;) Eo ) 


3 3 
T; m1 (exp ((Ey—E;)/kT,)—1)? 2 ) (; [—F'(Eo/kT;)] kT; 





Ill. 


Whereas in the preceding section we have 
discussed the properties of B.E. gases whose 
atoms were supposed to have arbitrary individual 
energy spectra, we should like to discuss now 
briefly the behavior of such gases in periodic 
potential fields. The internal structure of the 
atoms will be left aside. 

We should like to consider first the case of the 
so-called free-particle approximation. Here the 
amplitudes of the field’s potential are so small 
that the influence of the field on the motion of the 
particles consists essentially in the introduction 
of the well-known discontinuities in the energy 
spectrum of the individual atoms. We will hardly 
influence the behavior of our B.E. system by 
assuming that the free-particle level spectrum is 
realized everywhere in the interior of the per- 
mitted energy bands. The small but physically 
important deviations from this spectrum in the 
vicinity of the discontinuities will be neglected in 
order to simplify the calculations. We denote by 
e; the discontinuity in the energy of a particle 
when crossing the boundary between the ith and 
(t+1)th Brillouin zone in wave-vector space. We 
will have the foliowing distribution functions for 
our gas occupying a volume V at temperature T 
in the periodic potential field assumed to be 
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simple cubic with the spatial period d, 
V dpdp,dp: 


dni+1 =— 


h’ exp (a+fe:+B(p.?+p,?+p.7)/2m) =f 


and the parameter a determined through the total number WN of the gas atoms will be given by 


(29) 





N= 7 dnis1 
i=0 


V dpdp,dp. 
af 
h? i=o » €xp (a+Ber+B(p.?+ wii 


(30) 





-)/2m)—1 


where v is the fundamental volume in momentum space associated with each Brillouin zone. It is 
clear that €9 is zero, and that v is the volume h*/d*; the volume associated with the first zone being 
that of a cube centered at the origin of momentum space. In order to simplify further the calculations 
without reducing the generality of the discussion it will be assumed that there are but two zones, the 
first and the second, and from then on the particles may be regarded as being completely free. A 
simple calculation shows then that the condensation temperature of the system in this periodic field is 
the root 7; of the equation 


(31) 





4) F(a=0, €:/RkT;, €2/kT;, h/d(2mkT;)*) 
a F(0) 


where F(a, €:/kT, - ++) generalizes the often-met F functions associated with the free-particle continu- 
ous spectrum, and 7» is the transition temperature of the free gas with the same concentration in 
absence of the field. Since F(a, €:/kT, ---) < F(a) for all a, 7; must be necessarily larger than T>. 


The average energy of the gas will now be given by 


3 €2 h 
E_ == emer yere(o, * —,— i); 


<T;3; (32) 


kT kT d(2mkT)! 


3 V €1 
E,=- __(2rmkT)*TG( a — 
2 h3 kT 


where G(a, €:/kT, ---) generalizes the function 
G(x) associated with the free-particle continuous 
spectrum. The explicit expressions of these more 
general F and G functions are somewhat compli- 
cated and will not be given here. A simple glance 
on the forms of E_ and E, indicates that here 
too the order of the transition will depend on the 
value of (da/dT)r; since the heat capacities 
(dE_/dT)r; and (dE,/dT)r, differ only by a 
term proportional to (da/dT)7;. One can show 
now without difficulty that this derivative and, 
therefore, AC, vanish at the transition tempera- 
ture. The band structure of the individual energy 
spectrum with the free-particle level density does 
not change the nature of the transition from that 
existing with the completely free B.E. gas. This 


); T>T;; (33) 





result is due to the fact that around the origin of 
the momentum space, inside the volume corre- 
sponding to the first Brillouin zone, the atoms 
behave as free particles. It is clear that the 
approximation made in the expression of the 
energy and level density in the neighborhood of 
the energy discontinuities cannot affect quali- 
tatively this result. 

We should like finally to comment briefly on 
the B.E. transition in the limiting case of tight 
binding, i.e., of tightly bound B.E. particles in 
molecular periodic fields. The individual energy 
spectrum is formed now by a series of bands. One 
can show now that the condensation proceeds in 
this case also as in the more simple case of the 
free-particle approximation. However, the nature 
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of the transition will depend now on the level 
distribution near the bottom of the deepest 
energy band. In other words, the order of phase 
transition will depend now on the behavior of the 
bound B.E. particles around the origin of the 
first Brillouin zone. 


IV. 


The above given discussion of B.E. systems, 
apart from its academic interest has been 
motivated by the possibility, pointed out by 
London,” that the curious transition observed in 
liquid helium might be connected to the B.E. 
condensation phenomenon. On the basis of the 
preceding study it appears possible to make a 
few remarks on any model of liquid helium 
founded on a B.E. assembly which admits a 
distribution function. 

As we have seen, it seems possible eventually 
to define such conditions where the B.E. con- 
densation is accompanied by a discontinuity in 
the heat capacity. Neither the energy nor the 
entropy have anomalies in the cases we have 
considered. The anomaly in the heat capacity is 
connected to the fact that the derivative, with 
respect to the temperature, of the fundamental 
parameter a is finite at the transition point. 
Those quantities which do not depend on this 
derivative show no anomaly at the transition 
temperature. 

A model reproducing almost quantitatively 
several features of liquid helium, under its 
saturated vapor pressure, has been devised 
recently by London.* In order to get the ob- 
served discontinuity in the heat capacity it was 
assumed that the helium atoms, although having 
their energy proportional to the square of their 
momentum, had a level density which was 
proportional to a power of their momentum 
higher than two. Apart from the difficulty which 


3F. London, J. Phys. Chem. 43, 49 (1939). 
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might exist in any attempt to justify the associ- 
ation of such a level density to practically free- 
particle energies, there are two remarks which 
can be made about such models. 

First, the model in question has left out of 
account the physical condition which must be 
fulfilled by any B.E. gas model, namely, that for 
high energies the particles will necessarily behave 
practically as free particles, and at such energies 
the level density must become equal to that of 
free particles. 

The second remark concerns the preceding 
model as well as any other similar model. This 
can be formulated as a proposition which can be 
proved easily. The proposition in question is the 
following: Consider any B. E. system in which it 
is possible to define in a non-ambiguous way a 
distribution function. The condensation tempera- 
ture of such a system will always vary in the 
same direction as the total number of particles, 
the concentration, the density or the pressure of 
the system. The derivative, with respect to the 
temperature, of the functions representing these 
physical quantities are always positive along the 
transition curve for such B.E. systems. The proof 
of this proposition is quite simple and will be 
omitted here. 

This characteristic behavior of such B.E. 
systems contrasts strikingly with that observed 
in liquid helium where the mentioned physical 
quantities vary in opposite direction to the 
temperature along the transition or lambda line. 
One is therefore led to the conclusion that 
investigations of non-ideal B.E. systems are 
necessary. In such systems no individual level 
spectrum and no distribution function can be 
defined anymore. It seems that as long as the 
B.E. condensation in such non-ideal assemblies is 
not known, the question of the possibility of 
connection between the transition observed in 
liquid helium and the B.E. condensation phe- 
nomenon remains entirely open. 
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Methyl- and ethyl-induced polymerizations of ethylene and propylene were carried out in 
static, as well as flow experiments. The latter allowed collection of enough reaction products for 
analysis by low temperature distillation. Efficiency of polymerization in propylene is much 
lower than in ethylene which is explained by the easy formation of allyl radicals. Analyses of the 
reaction products readily indicate that the polymer is formed by the successive addition of 
olefin molecules to radicals, that chain transfer through radical decomposition takes place, and 
that the chains are terminated by disproportionation or association of two radicals. 





HE recent continued interest in free radical 
reactions! has prompted us to publish and 
discuss these data on the polymerization of 
ethylene and propylene by organic free radicals 
which were obtained several years ago and which 
appear to be of considerable interest in the light 
of recent publications in this field. 


APPARATUS AND EXPERIMENTAL METHODS 


Both static and flow experiments were carried 
out. For static experiments a bulb (100 cc) was 
immersed in a well-stirred molten salt bath of 
the required temperature, and pressure changes 
with respect to time were taken by means of a 
Bourdon gauge, the volume of which was small 
as compared to the volume of the bulb. To 
eliminate the effect of temperature changes on 
the gauge it was maintained at constant temper- 
ature in a glass jacket through which water was 
circulated. Ethylene, propylene, and all poly- 
merizing agents were freed of dissolved oxygen 
by solidifying or liquefying them with liquid air, 
distilling off a part of the material under vacuum, 
and using the center cut for the experiments. 

In the flow experiments the olefin was run 
through a flow meter, a drying tube, and into an 
evaporator containing either tetramethyl or 
tetraethyl lead. The temperature of the evapo- 
rator bath was so regulated that the desired 
quantity of lead alkyl was carried by the olefin 
into the preheater and furnace. The reactor was 
a Pyrex tube about 30 cm long. In preliminary 
runs the diameter was varied widely in order to 

1J. C. Jungers and L. M. Yeddanapalli, Trans. Faraday 


Soc. 36, 483 (1940); H. S. Taylor and J. O. Smith, Jr., 
J. Chem. Phys. 8, 543 (1940). 


detect possible effects of the wall but no influence 
was observed. The exit gases passed through 
tightly packed glass wool to remove lead; the 
polymers were condensed in traps cooled by a 
CO.-alcohol mixture; and in most cases the 
remaining gases were vented. 


EXPERIMENTAL RESULTS 


The static method was used mainly to investi- 
gate the relative efficiency of a great number of 
polymerizing agents. Of particular interest is the 
greater ease of the free-radical-induced polymer- 
ization of ethylene as compared with that of 
propylene.” Figure 1 shows the polymerization 
of ethylene and propylene by lead tetramethyl 
in which the decrease of pressure is plotted 
against time. The experiments were started 
slightly below atmospheric pressure. 
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2See also F. O. Rice and E. Teller, J. Chem. Phys. 6, 
489 (1938) who mention unpublished work by D. V.- 
Sickman. 
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Pb(C2H5)« 
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Polymerization reactions of this type do not 
lend themselves readily to a theoretical treat- 
ment on the basis of experimental reaction rates 
and yields because the different stages in the 
polymerization need the postulation of hypo- 
thetical steps which, preferably, should be 
verified experimentally. For this reason detailed 
knowledge of the reaction products appears to 
be of interest. Since relatively large quantities 
of reaction products were needed for analysis by 
Podbielniak distillation, the flow method de- 
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Fic. 2. Analysis of polymerization products (paraffins 
oe. For type of reaction see corresponding numbers 
in Table I. 


scribed above was adopted. The temperatures 
were chosen as low as was compatible with 
complete decomposition of the lead alkyls used. 
An analysis from an experiment using 4.06 mole 
percent lead tetramethyl showed 0.1 g of lead 
remaining in the product from 38.2 g lead alkyl 
consumed, showing the lead alkyl was practically 
completely decomposed. The conditions of this 
particular experiment were considered the most 
likely to show undecomposed lead compounds. 

Table I gives 6 representative experiments. 
The relative effectiveness of the lead compound 
in each experiment is denoted by the net moles 
of olefin combined per mole of radical consumed. 
It is seen that this yield drops off considerably 
for higher radical concentrations and that the 
yield is higher in polymerizations induced by 
methyl radicals than by ethyl radicals. As in 
the static experiments with methyl radicals the 
yield obtained with ethyl radicals in the poly- 
merization of propylene is low as compared with 
the ethylene polymerization. 

The analyses of the reaction products heavier 
than the gases to be polymerized are shown in 
Fig. 2. The fractions are plotted in percent of 
total polymerization product against the carbon 
number of the fractions. In Fig. 3 the analytical 
results of the methyl-induced ethylene polymer- 
ization (Exp. 1 and 2) are split up into paraffins 
and olefins. In Fig. 4 the same is done for the 
ethyl-induced ethylene polymerization (Exp. 4) 
and in Fig. 5 for the ethyl-induced propylene 
polymerization (Exp. 6). Some further segration 
is given in the figures. In one case—the ethyl- 
induced ethylene polymerization—a_ careful 
analysis of the gases not condensed in the CO, 
traps gave no hydrogen and only relatively small 
amounts of methane and ethane. 
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DISCUSSION 


The strikingly low polymerization efficiency 
of both methyl and ethyl radicals in propylene 
as compared to that in ethylene is at once 
explained in the light of recent experiments by 
H. S. Taylor and J. O. Smith.'! These investi- 
gators have measured the methane formation in 
the interaction of methyl radicals with propylene 
and find the surprisingly low activation energy 
of 3.1 kcal. for this process—the lowest activation 
energy found in any of their studies with 
hydrocarbons. The relatively very low yield of 
methane and ethane during the ethylene poly- 
merization in our experiments, together with 
earlier observations’ that a much_ smaller 
methane formation obtains from reactions of 
methyl radical with ethylene than with saturated 
hydrocarbons, suggests strongly that in our case 
the polymerizing radical chain (to be discussed 
below) in propylene may easily be broken by 
methane and ethane formation leaving an allyl 
radical.* The allyl radical having resonance 
characteristics and an estimated resonance 
energy of 15.4 kcal. would have a long life and 
would mainly disappear by collision with itself or 
other radicals, the latter having a net effect of 
breaking two radical chains. It is for these reasons 
that the interaction of methyl and ethyl radicals 


3H. S. Taylor and J. C. Jungers, Trans. Faraday Soc. 
33, 1353 (1937). 

* Recent mass spectrographic studies in these labora- 
tories by Mr. G. C. Eltenton on the interaction of alkyl 
radicals and olefins have shown the above formation of the 
allyl radical for the first time in a direct manner. This mass 
spectrographic work will be published in the near future. 
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with ethylene leads more easily to the formation 
of higher alkyl radicals. The analysis of the re- 
action products readily indicates that the polymer 
is formed by successive additions of the olefin 
molecule to radicals through a bimolecular re- 
action of the type 


CH3;+C2HiC3H:, 
C3;H;+C.H.>CsHu Pe ee 


The fact that most of the polymer is of low 
carbon number and that the products of low 
molecular weight contain more olefins than 
paraffins shows that the formation of larger 
radicals is mainly interrupted by decomposition 
of the radicals into olefins and smaller radicals, 
i.e., by chain transfer through radical decompo- 
sition. Other experimental evidence to this effect 
is found in a paper by Paneth and Lautsch‘ where 
the preparation of free propyl and free butyl 
were attempted by decomposition of tetra-n- 
propyl lead and tetra-isobutyl lead at 400°C 
and 300°C, respectively. Methyl radicals only 
were detected in both cases and the conclusion 
was drawn that free propyl and butyl were 
formed transiently, but very rapidly decomposed 
into ethylene and methyl or propylene and 
methyl. F. O. Rice and K. K. Rice in their 
monograph on free radicals’ also assume in 


iJ 
° 





Fic. 4. Analysis of 
products of ethyl-in- 
duced ethylene poly- 
merization (———— 
olefin, par- 

‘affin). percent 

C,Hs=0.3 percent iso- 

C4Hs+16.3 percent 

n-C4Hs+0.7 percent 

1,3-butadiene. 





a 





aw 


























° 


C4 Cs Ce 
CARBON NUMBER 


% WEIGHT OF TOTAL POLYMERIZATION PRODUCT 
3 


general that larger radicals break down very 
readily into smaller radicals and compounds 
containing a double bond. It is seen from 
experiments of Frey and Hepp* that about ? of 
the butyl radicals formed by decomposition of 
di-n-butyl mercury at 490°C decompose accord- 


4 Paneth and Lautsch, J. Chem. Soc., p. 380 (1935). 

5 F. O. Rice and K. K. Rice, The Aliphatic Free Radicals 
(The Johns Hopkins Press, 1935). 

6 Frey and Hepp, J. Am. Chem. Soc., 55, 3357, 1933; 
(see reference 5, p. 65). 
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ing to the equation: 


C,H g—C2H a+CoH 5. 


From the distribution of the products in the 
reaction between ethylene and ethyl radicals 
(Exp. 3) we must deduce that this and similar 
reactions leading to butylene or hexylene forma- 
tion occur simultaneously with reactions leading 
to propylene and amylene and CH; radicals. 

The paraffins must be formed either by 
disproportionation or by association of the 
radicals, i.e., by chain termination, for instance, 


2C4Hy—C4Hs+ CH 0 
2C,H 9—C;3H 18- 


The decomposition of an alkyl radical into a 
hydrogen atom and an olefin, which has been 
frequently suggested, does not take place since 
no trace of hydrogen has been found in our 
reaction products. 

Characteristic of the three types of polymer- 
ization performed is the fact that the efficiency 
increases with decreasing amounts of radicals. 
This is easily understood because there is less 
chance for two radicals to be eliminated by 
reaction with each other to form stable mole- 
cules, a process which furthermore is energeti- 
cally (dissipation of reaction energy) less favor- 
able for methyl than for ethyl radicals. This 


PROPYLENE 





r 
° 





a 


Fic. 5. Analysis of products 
of ethyl-induced propylene 
polymerization (——-—— ole- 
fin, paraffin). 8.9 
percent C,Hs=2.3 percent iso- 
C,Hs+6.6 percent n-C,4Hs. 








S 





wa 











° 





C4 Cs 
CARBON NUMBER 


% WEIGHT OF TOTAL POLYMERIZATION PRODUCT 


may account for the latter’s lower efficiency. 
Finally, it is of interest to note that in the 
reaction of ethyl radicals with propylene 35 
percent of the butylene formed is isobutylene. 
Unquestionably the only major reaction leading 
to its formation must be: 


C;H,+C2H;—iso C;H:;—iso C,H;+CH 3- 


This is in sharp contrast to the reaction of 
propylene with HBr for which a chain mechanism 
with bromine atom as carrier’ has been generally 
accepted and where the bromine atoms almost 
exclusively go to the end position. 

7M. S. Kharasch, H. Engelmann and F. R. Mayo, 


J. Org. Chem. 2, 288, 400, 577 (1937); F. R. Mayo and 
C. Walling, Chem. Rev., 27, 372 (1940). 
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Statistical Thermodynamics of Several 
Halomethanes 


WALTER F. EpGELL AND GEORGE GLOCKLER 
Department of Chemistry, State University of Iowa, Iowa City, Iowa 
April 16, 1941 


HE heat capacities of methyl fluoride and methyl 

iodide at several temperatures have been calculated 
to the rigid rotator simple harmonic oscillator approxi- 
mation in connection with another problem. Inasmuch as 
they do not appear elsewhere in the literature it was felt 
worth while to present them here. 

The fundamental frequencies are well known and appear 
in Table I together with those of CH2Bre2 as discussed 
below. The value 2916 cm™ is the mean of the resonating 
frequencies 2862, 2965 and 2860, 2971 appearing in methyl 
fluoride and methyl iodide, respectively. Linnett! has used 
the values 2862 and 2940 in calculating force constants for 
these molecules, but it appears more probable that the 
unperturbed fundamental lies closer to the mean. The heat 
capacities are given in Table II. Entropies at 298.1°K have 
also been calculated and appear in Table III. Moments of 
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inertia used were those given by Dennison.? The standard 
state is the usual one of unit fugacity and zero pressure. 
Corrections to finite pressures have been adequately 
discussed.?4# 

When the heat capacities of CH;Cl, CH2Cl, and CHC; ' 
(in the temperature range 7 = 298.1 to 600°K) are plotted 
against log 7, the resulting curves are nearly parallel. 
They may be brought into near superposition by shifting 
the temperature scale for each substance, i.e., plotting 


TABLE I. Fundamental frequencies. 








CH3F* CHs3I* CHoBre*-t 
2987 (2)** 3074 (2) 3060 
1476 (2) 1445 (2) 2988 
1200 (2) 885 (2) 1390 
2916 (1) 2916 (1) 1183 
1460 (1) 1252 (1) 1090 
1048 (1) 532 (1) 806 
636 
576 
174 











* H. Sponer, Molekiilspektren (J. Springer, Berlin, 1935). 
** Values in parenthesis are the degeneracies. 
t See reference 6 of text. 


TABLE II. C°p for several halomethanes. 








CH3F 


8.95* 

9.72 
10.09 
10.55 
11.80 
12.25 
13.83 











* Units are cal./°K/mole. 


TABLE III. S°298.1 for several halomethanes. 








CH3F CHsl CH2Brz 





53.30 E.U. 60.85 E.U. 70.16 E.U. 
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Fic. 2. Shift in temperature scale vs. molecular weight. A—CHsF, 
CH2F2, CHF3;; B—CH3Cl, CHeCle, CHCl3; C—CHs3Br, CHeBro, 
CHBr3; O—-CH2Brz, Stevenson and Beach.‘ 


C°, against log 7/6, where @ is a constant characteristic 
of each substance (see Fig. 1). If log 6 is plotted against 
molecular weight, the three points lie on a straight line. 
This is shown in Fig. 2. The same thing is true of the 
corresponding fluorine compounds. However the heat 
capacity of CH2Bre, as calculated by Stevenson and 
Beach*” cannot be brought into superposition with that of 
CH;Br and CHBr;—being too high at low temperatures 
and too low at high temperatures if the average value of @ 
is taken. This value of @ is designated by the clear circle 
in Fig. 2. Inasmuch as there is no question of the correct- 
ness of the fundamentals in these latter two molecules, it 
may be assumed that the fault lies with those of CH2Bro. 
Since the above calculations” were made, Wagner’s in- 
teresting normal coordinate treatment of the vibrations of 
the methylene halides has appeared. Although this work 
still leaves some question as to the fundamentals in this 
class of compounds, it indicates that there is no room for 
five low-lying fundamentals in CH2Bre as assumed by 
Stevenson and Beach. This is confirmed by calculations 
made in this laboratory. Therefore the most uncertain of 
these, 478 cm, is deleted, and 1183 taken as the ninth 
fundamental. This frequency has been observed both in 
Raman and absorption spectra® and is not satisfactorily 
accounted for otherwise. With the fundamentals as given 
in Table I the calculated heat capacity of CH:2Br2 is 
brought into superposition with those of CH;Br and CHBr; 
by the value of @ falling on the straight line C in Fig. 2. 
The heat capacity and entropy of methylene bromide 
appear in Tables IT and III. 

'J. W. Linnett, J. Chem. Phys. 8, 91 (1940). 

*D. M. Dennison, Rev. Mod. Phys. 12, 175 (1940). 

*P. Fugassi and C. Rudy, Jr., Ind. Eng. Chem. 30, 1029 (1938). 

‘ References to calculated heat capacities other than those from this 
note are from the following sources: (a) CHsCl, CH2Cle, CHCls: R. D. 
Vold, J. Am. Chem. Soc. 57, 1192 (1935). (b) CHsBr and CHBrs: D. P. 
Stevenson and J. Y. Beach, J. Chem. Phys. 6, 25 (1938). (c) CH2F:: 

- R. Leader, “Studies in molecular structure,” Ph.D. Thesis, Uni- 
versity of Minnesota (1940). (d) CHF3: George Glockler and W. F. 
Edgell, J. Chem. Phys. 9, 224 (1941). 

5 J. Wagner, Zeits. f. physik. Chemie B45, 69 (1939). 

® See references to spectra of this molecule tabulated by Ta-You Wu, 


Vibrational Spectra and Structure of Polyatomic Molecules (China 
Science Corp., Shanghai, 1939), p. 248. 


485 


Thermodynamic Functions for Molecules with 
Internal Rotation 


KENNETH S. PITZER AND WILLIAM D. GwINN 
Department of Chemistry, University of California, Berkeley, California 
April 12, 1941 


EVERAL years ago one of us! presented formulas and 

tables of thermodynamic functions for molecules with 
restricted internal rotation. The cosine function was as- 
sumed for the potential barrier. The treatment was based 
on Nielsen’s* work on the quantum mechanics of two co- 
axial tops, and while not absolutely exact, it was a very 
good approximation for that case. It was presumed, on the 
basis of the classical approximation, that the results would 
probably be quite accurate for all molecules whose mo- 
ments of inertia are independent of the angle of internal 
rotation. A molecule of this class may also be described as 
consisting of a rigid frame with an attached symmetrical 
top. We can now prove rigorously this more general ap- 
plicability of the tables. In view of this development, we 
have decided to recalculate the tables with greater nu- 
merical accuracy, and to publish the new derivation with 
the new tables. The purpose of this communication is first 
to state the more general applicability of the old tables, 
and second to outline another, slightly less precise, deriva- 
tion which is based on the detailed analysis of the quantum 
mechanics of such molecules presented recently by 
Crawford.’ 

Crawford’s final formula is somewhat different from that 
used in computing the tables referred to above.* However, 
he introduced certain approximations which we shall 
attempt to remove. These approximations arose through 
the use, in correction terms, of quantities (v and a in the 
section “limiting values”) calculated with harmonic 
oscillator or rigid rotation wave functions, instead of the 
exact Mathieu functions or linear combinations thereof. 
While it is not possible to obtain a simple final formula us- 
ing exact correcting terms, it can be shown! that this final 
formula will depend only on the following independent 
variables: T, the potential barrier, V, the reduced moment 
of inertia F or J,, the number of potential peaks, v or n, 
and the ratio of the reduced moment to top moment {. 
If we can find the true partition function for some special 
case which includes all values of these variables we will 
have the problem solved simply and accurately. This 
is done by taking the case of two coaxial tops. Here 
Q=1,/(1;+J2) so that all values from 0 to 1 can be ob- 
tained. Likewise all values of the other variables are 
attainable in this special case. Since the tables! previously 
published were quite accurate for the case of two coaxial 
tops, we conclude that they are also quite accurate for the 
more general case of a rigid frame with an attached top. 
In the more general case Crawford’s formula for the reduced 
moment of inertia should be used. 


3 
= Foam al -2 asns| 


i=1 


Here the J;'s are the three principal moments of inertia of 
the whole molecule, the \;’s are direction cosines between 
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the top axis and these principal axes, and A is the moment 
of inertia of the top itself. 

To summarize briefly, Crawford’s analysis shows that 
the partition function does not change when one goes from 
the special case of two coaxial tops to the more general case 
of any rigid frame with one top. Therefore one may use the 
simple tables computed for the special case in the more 
general case. The problem of a frame with more than one 
attached top will be discussed in the paper to be published 
later. 

This result affects the conclusions of a number of recent 
publications. The general recommendation of Wilson® that 
Crawford’s formulas are more accurate is to be reversed. 
Insofar as the two methods differ, the new calculations of 
Crawford, Kistiakowsky, Rice, Wells and Wilson’? on 
propylene are in error. Since we have believed that their 
conclusions were uncertain on other grounds,® we merely 
conclude that doubt remains. In Crawford’s® calculations 
for methyl alcohol, the results based on his formulas should 
be discarded in favor of results from the tables which he 
also gives but dismisses as less accurate. This reduces, but 
does not remove, the discrepancy between the entropy 
calculated from the 470+40 cm™ barrier proposed by 
Koehler and Dennison!’ and the experimental value. We 
may also note that Aston and Doty’s" calculations on 
methyl amine may now be said to be based on accurate 
equations. Consequently one concludes that the torsional 
frequency in methyl amine is greater than the 270 cm™ 
value suggested by Owens and Barker.” 

1K. S. Pitzer, J. Chem. Phys. 5, 469 (1937). 

2H. H. Nielsen, Phys. Rev. 40, 445 (1932). 

3B. L. Crawford, Jr., J. Chem. Phys. 8, 273 (1940). 

4 For the extensively investigated case of ethane, Crawford's formula 
is a much poorer approximation that that used in computing the tables. 

5 This functional dependence can be verified easily for Crawford's final 
approximate formula, Eqs. (2-7), or for his accurate formula through 
the first order correction, Eqs. (41) and (42). To obtain this result accu- 
rately through the second-order correction, it is necessary to note that 
his quantity (H2)* ¢ will still be proportional to pm? but with a coefficient 
in general different from —Fm/2. This coefficient, however, will be a 
function of Vm, vm, Fm, and ¢ only. One must carry out a transformation 
similar to that between the unnumbered equations at the top of p. 280, 
column 2. The desired result may then be verified without great 
difficulty. 

6 E. B. Wilson, Chem. Rev. 27, 17 (1940). 

7B. L. Crawford, Jr., G. B. Kiistiakowsky, W. W. Rice, A. J. Wells 
and E. B. Wilson, Jr., J. Am. Chem. Soc. 61, 2980 (1939); also J. 
Chem. Phys. 8, 610, 970 (1940). 

8 K.S. Pitzer, Chem. Rev. 27, 39 (1940). 

9 B. L. Crawford, Jr., J. Chem. Phys. 8, 744 (1940). 

10 J. S. Koehler and D. M. Dennison, Phys. Rev. 57, 1006 (1940). 
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The Photovoltaic Effect 


EpWARD ADLER 
Department of Chemistry, College of the City of New York, 
New York, New York 
April 12, 1941 


HIS is in reply to a letter by H. M. Cassel about our 
papers! on this subject which appeared in the April, 
1941 issue of this journal. 


THE EDITOR 


We do not wish to enter into a polemic on this subject, 
but we believe there are some statements in the letter 
which are not entirely valid, and which require an answer. 

(1) The theory which we proposed for the photovoltaic 
or Becquerel effect enables one to derive a simple mathe- 
matical relation between light intensity and potential 
which is in agreement with the experimental light in- 
tensity/potential curve for small and medium light in- 
tensities. The term accounting for the photo-current, due 
to the transition of photoelectrons from the electrode to 
the electrolyte: 

i, =ek;J (1) 


(using our previous nomenclature) is mot introduced ad hoc 
into an overvoltage equation, as Dr. Cassel’s letter seems 
to imply, but follows naturally from the concept of a 
potential barrier. 

The fact that our equations are similar to those en- 
countered in the treatment of overvoltage means simpl\ 
that similar processes are present in both phenomena, and 
that both problems can be treated by the method of 
absolute reaction rates. 

(2) Our derived relation referred to above does not, 
indeed, without modification lead to the ‘saturation po- 
tential’’ at very high light intensities, which is experimen- 
tally found. We were, of course, well aware of this so-called 
“‘discrepancy”’ and suggested that it might be due to the 
fact that only a finite, limited number of photo-active 
centers were available for liberating photoelectrons. 

Dr. Cassel believes that a ‘‘saturation potential’’ may be 
explained by the substitution for (1) of 


I, =€[Atads Jhs’J, (2) 


and indicates that “It can easily be seen that with this 
correction the theory covers the whole range of the em- 
pirical curve.” We do not believe that this can be at all 
easily seen, for if the substitution zs actually made and the 
mathematical operations carried through, the result 7s 
exactly the same as the one which we found for medium 
light intensities, namely, Ve= C+B log J. 

(3) The fact is, of course, that the number of photo- 
active centers is finite and limited. This is not an assump- 
tion. It can be seen from data in the field of photo- 
conductivity, and also in dry, solid photo-cells in both of 
which a saturation current and a saturation photo-potential 
obtain. It seems to us that in the photovoltaic effect also, 
at very high light intensities, the limited number of photo- 
centers would make itself noticed in the form of a satura- 
tion effect. We might also indicate that in our paper we 
pointed out that a primary negative photo-e.m.f. could be 
explained by a mechanism into which [Atads] does nol 
enter and which gives the same type of relation found for 
primary positive photo-e.m_f. 

1 E. Adler, “The photovoltaic effect,’’ J. Chem. Phys. 8, 500 (1940); 


“The photovoltaic effect,’’ Trans. Electrochem. Soc. Preprint 79, 6 
(1940). 
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